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Introduction and plan of the notes 



These notes have been written for a short course at the VIII Modave school of Math- 
ematical Physics, and aim to illustrate some aspects of universality and Renormal- 
ization Group (RG) techniques in Physics. 

The topic is broad and complex. Even focusing on a particular aspect would 
take a whole monograph, and indeed there are several dedicated to this purpose. 
Here the focus is on the relatively new approach of Functional Renormalization 
Group Equations (FRGE, sometimes called, misleadingly in our opinion, "exact"), 
and in particular on Wetterich's equation, which is being widely applied in statistical 
Physics, Standard Model Physics, and to study the asymptotic safety conjecture of 
Weinberg for Quantum Gravity. 

The course from which these notes originated was intended for young re- 
searchers that are not necessarily familiar with renormalization group ideas, a topic 
that often finds little space in master courses in theoretical Physics. Therefore, 
rather that focusing on Quantum Field Theory (QFT) alone and going deep into 
the many technical tools needed to employ the FRGE in the current topics of re- 
search, we preferred to keep the discussion as elementary as possible. In this spirit, 
the ideas of universality and renormalization group are first presented in the simple 
context of dynamical systems. 

The first part of the notes introduces the basic notions related to dynami- 
cal systems, including fixed points, attractors and asymptotic manifolds, which are 
crucial for understanding FRGE in quantum field theory. Then Feigenbaum's uni- 
versality for iterated maps on the interval is presented, and explained in terms of 
renormalization group, following the exposition of Collet and Eckmann. This has the 
advantage of introducing RG in a relatively simple context, where the mechanism 
has been completely understood and mathematically formalized. 

In the second part, we focus on QFT, restricting to the case of scalar fields. 
First we briefly recall some textbook notions about QFT: divergences, renormal- 
ization, anomalous scaling, etc., flxing the notation for the following sections. The 
we present the formal derivation of Wetterich's FRGE and apply it to obtain again 
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the perturbative /3-functions of theory, commenting on the differences with (di- 
mensional regularization) perturbative calculations. Finally, we use the FRGE to 
investigate the Wilson-Fisher fixed point for three-dimensional Ising-like magnets, 
a calculation quite reminiscent of the one for the Feigenbaum fixed point. 

Many things are left out of this presentation. The general spirit is to present 
all ideas in the simplest possible context, and extend them to more realistic scenarios 
by analogy. The mathematical formalization is minimal by necessity of elucidating 
the more conceptual points. This can be a bit unsatisfactory in particular in the 
first part, where the such formalization could be used to provide rigorous proof of 
universality. In the second part, many important topics are omitted, such as the ap- 
plication of FRGE to theories with (gauge) symmetries and on curved backgrounds, 
as well as a discussion of the asymptotic safety conjecture for gravity; however, such 
topics may be found elsewhere in the literature. Good companions of these note 
are the ones of two parallel Modave courses that introduce conformal field theories 
[63 ] and elucidate the relation between RG flows and scale invariance in QFT 
topics that unfortunately do not find space here. 

Finally, here the FRGE is presented as an useful tool but no statements are 
made about its reliability for rigorously investigating non-perturbative scenarios and 
establishing renormalizability in a constructive sense, which we believe is yet to be 
demonstrated. We hope to return to these issues in the future. 

To make the presentation more pedagogical, several figures are included, as 
well as a selection of short exercises. 
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Part 1 

Dynamical systems 



Dynamical systems constitute a branch of mathematical Physics concerned, roughly 
speaking, with study of the time evolution of some space. They can be used to 
study many types of problems, from chemical reactions to fluctuations of the stock 
market, and including all of Hamiltonian mechanics. 

The type of time evolution (deterministic or stochastic) and the nature of the 
space under consideration (a smooth manifold, a measure space, etc.) dictate the 
kind of tools to be employed in this study. Here we will consider only the simplest 
type of dynamical systems, where evolution is deterministic and the space is some 
(open subset of) M". General references references are e.g. [TOl [211 121 162] . 

These simple cases can be split in two categories depending on which notion 
of time we decide to adopt that is, whether the system will have continuous time 
t G M or discrete time t G Z. Let us start from the former. 

1.1 Dynamical systems with continuous time 
1.1.1 Generalities 

In this case we are with dealing autonomous, flrst order ordinary differential equa- 
tions (ODEs) on an open set Ai G M", of the form 

(1.1) x = f{x), xeM. 

To keep our discussion as simple as possible, we assume everything to be suitably 
regular, so that Cauchy's theorem guarantees the existence of the unique solution 
with initial condition x G at time t = 0. We will call this solution $*(a;). 
Furthermore we will assume that this solution exists for all t G M (or at least for 
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all t > 0). In fact, in what follows, we will not be interested in solving one specific 
Cauchy problem, but on understanding the generic motion of a generic point x G A^, 
and in particular in what happens asymptotically, i.e. when t — )■ oo. To this end, 
we can study the map 



(1.2) M^M 



which will be smooth under our assumptions. In the case where the solution exists 
for all t G M, one immediately notices the following properties 



(1.3) <I)° = Id, ($*)-! = $-*, <|.*o$^ = $*+% 

that imply that we can define a one-parameter abelian group of diffeomorphisms 

(1.4) $ = {$*, tGM}. 

One could also distinguish the case where $(a;) is not invertible, and consequently 
$ is only a semigroup, but this is not important now. We will say that $ is the flow 



of the differential equation (1.1). 



We are now in a position to give a more formal definition. A continuous time 
regular dynamical system is a couple {Ai, $) where Ai i^a regular open subset of 
M" and $ is a one parameter group of diffeomorphisms. 

Before proceeding further, it is worth pausing for some comments. We said that 
we will be mostly interested in the asymptotic behavior of these systems. Why? Let 
us be concrete, and pick the study of the solar system as a specific case. In principle, 
we could solve Cauchy's problem for given initial conditions on the positions and 
velocities of the planets. In practice this is impossible analytically, but for the sake 
of the argument let us assume that we can do it somehow (e.g. numerically) to any 
desired accuracy. Then, since Newton's equations are deterministic, we should know 
exactly the orbits of any planet at any instant]^ 

In modern times, after the works of Poincare, this idealization has been con- 
sidered less and less convincing, for reasons that are probably clear to the reader. 
First, we cannot solve any ODE to any accuracy, and besides, when we decide to 
treat the solar system as isolated we overlook the possibility of a small perturbation 
occurring at some point (e.g. a small asteroid passing through it) that would alter 
the equations. Even if we could eliminate these problems, we still have to specify 
some initial conditions, that come from a physical measurement, and carry some 



^More generally, Ai can be a n-dimensional smooth manifold. 

^ This reasoning could be called Laplacian determinism, as the XIX century mathematician 
Laplace was among the ones that emphasized it most. 
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uncertainty, i.e. x{t = 0) = Xq ± Sx. Using regularity of the ODEs, one can find a 
bound on the propagation of this uncertainty: schematically 

(1.5) \\^'{xo + 6x) -^\xo)\\<Ce^'6x, C>0, A>0. 

In practice, due to the wild exponential behavior, even when 6x is very small, the 
separation will become large after a short interval of time. 

It is for these reasons that it is more interesting to study the properties of 
$(x) for a generic a; G A^, and ask asymptotic and often qualitative questions, 
such as whether the solar system is stable, or some planet will be ejected from it at 
some point. This is an important point to keep in mind, not only in the context of 
dynamical systems, but of Quantum Field Theories (QFTs) as well. 

Exercise 1.1 (First order ODEs). We are restricting here to (systems of) first 
order, autonomous ODEs, which could seem to be a strong limitation. On the con- 
trary, show that an (autonomous or not) ODE of any order can be recast as system 
of first order autonomous ODEs. 

Exercise 1.2 (ODEs, solutions, and existence for all times). We are restrict- 
ing to ODEs such that $* exists and is unique for all t G M. One should realize that 
this is truly a restriction. 

1. Recalling Cauchy's theorem, write down a Cauchy problem that has infinitely 
many solutions. 

2. Consider the ODE with = M and given by x = 1 + x'^ . Does the solution 
exist for allt eR? 

Exercise 1.3 (Dependance on initial conditions). We want to illustrate the mag- 



nitude of the effects in (1.5). A simple, perhaps a bit trivial example, is a pencil 
balanced on its tip (inverse pendulum). The pencil has length L, and it is balanced 
on its tip with no speed and up to a great accuracy in the position of its top, 6. 
After how long a time At we won't be able to predict the position of the pencil with 
a precision of at least 1 rad ~ 57° ? 

You can take L = 10 cm, and try 6 = 10~'^m, S = 10~^^m. 
(Answer: At ^ 0.7s and At ^ 3.2 s.) 

1.1.2 Asymptotic behavior: fixed points 

Here we will discuss the asymptotic behavior for autonomous ODEs. The simplest 
case is when there exists one point x* E Ai such that 



(1.6) 



<!>\x*)=x*, Wt>T, 
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for some T. It is clear that this is equivalent to requiring the above condition to 



hold for all t and, in terms of (1.1), to requiring that 



(1.7) /(x*)=0. 

We will call such fixed point or critical point for the dynamical system. 

Of course it is not typical that the initial condition for a Cauchy problem is 
precisely x(0) = x*. These points are interesting because they influence the flow for 
any initial condition close to them. Consider, for instance, the case where = M 
and there exists an unique critical point x* such that f'{x*) ^ 0. Then, two things 
may happen: 



• f'{x*) > 0: then /(x) < to the left of x*, and the flow pushes these points 
to smaller values, away from x*; similarly f{x) > to the right of x* and the 
flow pushes them to the right, again away from x*. 

• f'{x*) < 0: then f{x) > to the left of x*, and the flow pushes these points 
to larger values, towards x*; similarly f{x) < to the right of x* and the flow 
pushes them to the left, again towards x*. Intuitively]^ the flow cannot cross 
X*, so that the motion tends to x*. 



We leave it to the reader that is new to these ideas to familiarize himself with this 
idea, e.g. by drawing more general scenarios for f{x) when Ai = 'R. 

This simple example motivates the need to classify fixed points depending on 
their property to attract or repel the points in their neighborhood under the flow 
$. This classiflcation, together with many useful criteria, was flrst put forward by 
Lyapunov. Let x* be a critical point for $, and let $ exists for all t G M. Then 



1. X* is attractive (or asymptotically stable) if there exists a neighborhood V of 
X* such that 

X eV =^ lim $*(x) = X* . 

t-5- + 00 

2. X* is stable for all time^if for any neighborhood U of x* there exists a neigh- 
borhood Vo of X* such that 

xeVo ^\x)eU VtGM. 



3. X* is unstable if it is not stable. 



^See Exercise 1.5 

^Stability only m the future or past amounts to restricting to t > t^^ or t < respectively. 
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4. X* is repulsive if there exists a neighborhood V of x* such that 

X e V ^ hm = X* . 

t— oo 

With this terminology, it is easy to classify the dynamical systems given by a 
linear ODEs. As it is well known, we have in that case 

(1.8) x = Ax, $*(x) = e*^a;, 

where e*"^ is the exponential of a matrix, defined by the convergent series e*"^ = 
Yllk=Q^^ / ■ Furthermore, if A is symmetric]^ then it can be diagonalized, and it 
is enough to consider its eigenvalues Ai, . . . , A„. Clearly, there is just one fixed point 
X* = 0, and it is not hard to verify the following statements: 

• If Re(Aj) < for all eigenvalues, then x* is attractive. 

• If Re(Ai) > for all eigenvalues, then x* is repulsive. 

• If Re(Aj) = for all eigenvalues, then x* is stable. 

• If Re(Aj) < for some eigenvalues and Re(Aj) > for others, then x* is 
unstable. 

We will also say that x* is hyperbolic if Re(Ai) ^ for all i. 

Of course linear equations are not very interesting per se. However, consider 
any ODE with (at least) one fixed point x*. Then we can write 

(1.9) X = fix) = A{x-x*) + 0{\\x- x*f) , with v4 = Jac/|^^^. , 

simply by expanding around x*. Our intuition suggests then that, at least if x* 
is hyperbolic, then the linearized analysis should be enough to "understand" the 
flow, at least close to x*. This is actually the it follows from an important 

theorem that we will state without proof (see e.g. [62]). 

Grohman-Hartman theorem: if i; = f{x) has a hyperbolic fixed point x* G 
Ai = MJ^, then there is a neighborhood of x* such that the flow there is homeomor- 
phic to the flow of the linear system x = Ax. In other words, locally the nonlinear 
flow is conjugated by a continuous invertible map to the linear one. There are ex- 
tensions of this theorem that guarantee (under additional hypotheses) that the two 
flows are diffeomorphic, but they are more subtle and we will not discuss them here. 



^More generally, similar considerations can be made using the Jordan form of A, but they will 
not be important for us. 
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Figure 1.1: Flow for a nonlinear and linearized system. 

As a corollary of this theorem, it follows that the stability property of an 
hyperbolic fixed point can be found from the corresponding linearized ODE, i.e. 
from the eigenvalues of the Jacobian at the fixed point, which is also known as 
Lyapunov's spectral method. 

Before proceeding further, it is worth illustrating all this on a simple example. 
Consider the following dynamical system, which we spell out in coordinates: 

(1.10) ii = a;i - xi a;2 , ±2 = -X2 + xj. 
One easily finds three fixed points 

(1.11) a;^i) = (0,0), x^2) = (l,l), x^3) = (-l,l)- 

Consider for instance x*^^y It is clearly hyperbolic and unstable, because the Jacobian 
there is the matrix A = diag(l,— 1). Therefore, by Grobman-Hartman theorem, it 
follows that the flow in the vicinity of x^^^ is conjugated to the one of the associated 
linear system. The latter is very simple: in the linear system, the axis xi supports 
an exponentially repulsive motion, whereas the axis X2 supports an exponentially 
attractive one. Generic initial conditions yield hyperbolae that asymptote to the 
coordinate axes. 



Looking at Figure 



1.1 



it is clear that the nonlinear flow around similar 



to the linear one. It is interesting to look for some curves that play a role similar to 
the coordinate axes in the linear system, which helped us to understand the motion 
of a generic initial condition. A very important and useful result guarantees their 
existence. 
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Stable manifold theorem (Perron-Hadamard) : ii x* is a hyperbolic fixed point, 
then the two sets 

(1.12) W = IxeM: \im(!>\x) = x* 

W = IxeM: hm (!>\x) = x' 

I t— >— oo 

are regular manifolds called stable and unstable manifold, and are tangent to the 
hyperplanes generated by the eigenvectors of the linearized system corresponding 
resp. to negative and positive eigenvalues (see e.g. [21 [38]). 

We shall not prove this theorem, but it is worth noticing that the difficulty 
is in proving that a regular stable (unstable) manifold exists locally. Once one has 
constructed such W^^^ (^loc) easy to obtain the whole manifolds by defining 

(1.13) = u ^\wu , = u $*(wroj . 

t<0 t>0 

We conclude this section with a word of warning: we discussed only the dy- 
namics around fixed points, but these are not the only objects that infiuence the 
asymptotic behavior. From dimension n > 2, dynamical systems may present limit 
cycles, and for n > 3 they may have chaotic behavior. Even giving an appropriate 
definition of these concepts will take us too far, and we invite the interested reader 
to consult, e.g. [70] . 



Exercise 1.4 (Qualitative analysis of ODEs). If the reader is not familiar with 
the qualitative analysis of ODEs, we suggest some exercises here. 

1. Consider the dynamical system x = Ax on M. = , where A is a symmetric, 
constant matrix, and draw the trajectories (the phase portrait) of the system. 

2. Consider a mechanical system of the form x = —V'{x). What are the stable 
and unstable fixed points in terms of V{x) ? Can the system have attractive 
fixed points? 

3. Perform the qualitative study of the Lotka- Volterra system 

(1.14) xi = a xi — (3 X1X2 , X2 = 7 2^2 + 5 a^ia;2 , (^,13,1,^ >0. 
You will find a detailed and commented analysis of this system e.g. in ITUV. 



Exercise 1.5 (Some short proofs). In stating our propositions we accepted some 
facts that the interested reader could prove, for a regular continuous time 

dynamical system. 
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Figure 1.2: Bifurcations diagrams showing the fixed points as functions of /x. 

1. Prove that if x,y G Ai, x ^ y, then ^ (t>^{y)- Conclude that a trajectory 
can never go past a fixed point. 

2. Prove that hyperbolic fixed points are always isolated. 

3. List reasonable properties that W^^^ and W^^^ should have, and assume them 
to prove the stable manifold theorem in two dimensions. 

1.1.3 Bifurcations 

We introduce here the notion of family of dynamical systems, that we will ana- 
lyze in more detail for discrete-time systems. Here we simply allow for a (regular) 



dependence of (1.1) on one or more real parameters yUi, . . . so that we have 



(1.15) X = f{x; fii, . . . ,fin) , xeM. 

This generahzation is quite natural, as it allows to study the evolution of a system 
depending on some external condition, such as the demographics as a function of 
resources, etc. 

Probably the simplest example is the linear equation 

(1.16) x = ijx, xeM = R, /xGM, 

and it is clear that the asymptotic properties of the system depend on the sign of /i. 
A similar, but somewhat less trivial example is e.g. 

(1.17) X = {x - - ij,2 , xeM = R, fii^R. 



1.2. Dynamical systems with discrete time — 9 



Here the parameter /xi is largely inessential, as it can be reabsorbed by a translation. 
However, again the sign of fi2 is important: when it is negative there are no critical 
points, whereas when it turns positive a couple of critical points (one stable, one 
unstable) is created. This is called the tangent bifurcation. It is useful to plot the 
bifurcation diagram (Figure 1.2(a)), where one draws the position of the critical 
points as a function of the relevant parameter (in this case, //2)- This gives an 
overview of the asymptotic properties of the system as the function of the external 
parameters. 

Another example, on which we will return later, is the so-called pitchfork 
bifurcation, given for instance by 

(1.18) x = -x^ + fix, xeM=R, /iGM. 

Here we have a stable critical point for < which splits up into two stable and 
one unstable critical points when n becomes positive. The bifurcation diagram of 
Figure 1.2(b) suggests the name. 

Exercise 1.6 (A bifurcation diagram). Draw the bifurcation diagram for 

(1.19) X = fi + x -]-x^, xeM = R, neR. 

Argue that, by changing fj,, this system can be used as a model of hysteresis in a 
magnet. 

Exercise 1.7 (Redundant parameters). Consider the logistic equatiorj^ 

(1.20) X = X - ij,2 x^ , xeM=R, /ii>0. 

Argue that by appropriate reseating, both parameters can be eliminated to yield 

(1.21) x = x-x'^, xeM = R. 



1.2 Dynamical systems with discrete time 

Some phenomena where observables can be measured only at some moment in time 
are more naturally described by using a discrete time variable. Examples are the 
abundance of a certain species after each reproductive cycle, the amount of crops 
collected every year, etc. 

To discuss dynamical systems with discrete time we just have to rephrase what 
we said in the previous section. 

^This describes the growth of a population with limited amount of resources, as opposed to 
Malthusian, or exponential, growth (see next section). 
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1.2.1 Generalities, fixed points, Lyapunov exponents 



Rather than being defined by an ODE such as (1.1), the typical definition of a 
discrete-time dynamical system is a recursion law 

(1.22) Xn+l = f{Xn) , f:M^M. 

We immediately obtain that the (discrete-time) flow $ satisfies, for n G N, 

(1.23) $° = Id, ^\x) = f{x), $-= ($i)" = /o...o/, 

where the last equation indicates the n-fold composition of functions. Depending 
on whether f[x) is invertible, one can add the additional property 

(1.24) $-1 = ($1)-^ = /-!, 



and extend (1.23) to n G Z, in which case the fiow will be a group, rather than just 
a semigroup. 

As seen in the previous section, it is interesting to look at fixed points x* E Ai 
that satisfy 

(1.25) <l''^(x*) = x* ^ x* = f{x*). 

Again, it will be important to understand whether x* attracts or repels the nearby 
points. Let us consider the case where = J C M is a (possibly unbounded) 
interval on the real line, and study the fiow oi Xq = x* + e. Then 

(1.26) xi = /(xo) ^x*+e fix*) , x„ = /" (xq) ^ x* + e (/'(x*))" . 

Clearly, the asymptotic behavior around x* depends on whether the modulus of the 
slope |/'(x*)| is larger than one (repulsive fixed point) or smaller than one (attractive 
fixed point). Similar notions can be formulated also when Ai is more general, i.e. 
in higher dimensions, but we will not need them in what follows; in fact, from now 
on, we will restrict to the simple case where Ai is an interval. 

Let us introduce a notion that will be useful later on, the one of Lyapunov 
exponent. This is a useful tool to understand the behavior of two neighboring [e- 
close) generic points in Ai: will they remain close together, and eventually get 
squeezed to the same attractive fixed point (as most of the cases that we have 
encountered studied so far), or will they become more and more separatee^? The 
natural quantity to consider is the the limit 

(1.27) 5(x,n) = lim^^^^^^t^ip^!M = |/'($n-i(a;))| . 



^ The latter is a typical feature of chaotic systems, which we will not define in detail. The 
idea is that in these systems the unpredictable behavior seen in Exercise [T3] for a particular initial 
condition (a pencil balanced on its tip) will be true for any initial condition. 
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To remove the dependence on n one can take the average of 5{x, n) along the orbit. 



Finally, since (1.26) suggests that the separation grows geometrically, so that we 
write 

N-l 

(1.28) 7W= lim -5^1og|/'(<l>"(x))| . 

n=0 

We will say (if the above limit exists) that 7(0;) is the Lyapunov exponent of x. A 
theorem by Oseledec [56j guarantees that indeed the limit exists for almost every 
X & Ai, and it is immediate to see that 7(0;) will be the same for any s G with 
the same asymptotic behavior. In the cases of our interest, in fact, there will be 
only one Lyapunov exponent, so that we will from now on drop the dependence on 

X. 

We conclude this section with two simple examples of one-dimensional discrete 
time-systems, that are also called iterated maps of the interval. As we already seen 
it is interesting to allow for dependence on one or more parameters /Xj. The simplest 
example is the Malthusian growth, a simple model for population expansion with 
unlimited resources. The law is linear 

(1.29) Xn+i = f^,{xn) , ff,{x)=fix, xe[0,+oo), fl>0, 
and the recursion can be solved immediately to give 

(1.30) $^(a;)=/x"x, neZ. 

As indicated, the flow can be extended to negative n, and dictates a simple geometric 
behavior, similar to the one of the continuous-time system x = log fi x. 

One could have the impression that discrete time systems are just a trivial 
modification of ODEs, but this is not the case. In fact, discrete-time systems given by 
a simple law can yield an extremely complex behavior. The model immediately more 
complicated than Malthusian growth is the logistic growth, which is a modification 
of the former to include finite resources. It can be written as 

(1.31) X n+i = f^,{xn) , f^{x) = fix {1 - x) , X G [0, 1] , 0</i<4, 

and a partial study of it, in its simplest regime, will occupy the rest of this chapter. 

This surprising complexity comes from the fact that discrete-time systems can 
be seen as arising from continuos-time ones in higher dimensions, by defining $"(a;) 
as the intersection of $*(x) with some submanifold embedded in A4, a procedure 
called Poincare section. It is then clear that discrete-time systems do not suffer 
the same topological limitation of continuous time ones, and therefore can exhibit a 
richer behavior also in low dimensioiS 

^To better understand what "topological limitations" we are referring to, the interested reader 



is invited to do Exercise 1.10 or read the proof of Poincare-Bendixon theorem in e.g. [70] . 
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Exercise 1.8 (Lyapunov exponent). Consider the iterated map on the interval 
M = [0, 1] defined by 

(1.32) Xn+i = f{xn) , fix) = 6{x - x^) (1 - 2x + 2x^) . 

What is the Lyapunov exponent? Check it numerically too. 
(Answer: log(2^/3 - 2).) 

Exercise 1.9 (Redundant parameters for the logistic map). The most gen- 
eral form of the discrete-time logistic map is 

(1.33) f^^^^^{x) = fiix - fi2x'^ , a; G [0,/ii/4/i2] , /ij > . 



Show that this can be reduced, by a rescaling, to (1.31), but the remaining parameter 
fi cannot be eliminated. Explain what is the difference between this case and the one 
of Exercise | i. 7[ 

Exercise 1.10 (Limit sets in D = 1). The set L^^ G Ai is called the limit set 
(more properly the cj-limit set) for $ if 



;i.34) L^ = \y: y= lim <l>*"(a;)| 



for some sequences {tnjneN and points x E Ai. Prove that is ^ is a continuos-time, 
one-dimensional dynamical system, then if L^^ ^ ^ it always contains a fixed point. 
Show with an example that this is not the case for discrete-time systems. 



1.2.2 The logistic map 

Let us start the study of the logistic map. As we said is defined by the one-parameter 
family of functions on = [0,1] 

(1.35) : [0, 1] [0, 1] , f,: x^^ix{l-x). 

In order to have that /^([0, 1]) C [0,1], it must be < ;U < 4. Let us start our 
analysis from small /i. 

When /i = 0, clearly nothing happens. When < /i < 1 the fixed point 
equation f^{x) = x has only one solution in x* = 0, which is a stable fixed point. 
For larger values of fi, x* = becomes unstable and a new stable fixed point is 
generated at x* = Its stability can be checked by looking at 

(1.36) |/'(x*)l = |2-/i|, 
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(a) /i — 0.9, equilibrium at x* — 0. (b) /i = 2.85, equilibrium at x* ^ 0. 




(e) n — 3.9, chaotic behavior. 



Figure 1.3: Asymptotic behavior of the logistic map for some /i's. 



so that this fixed point is stable for 1 < /i < 3. Intuitively, it is not easy to 
understand what happens when /x gets larger than 3; a way to proceed is to simulate 
the behavior of this system with a computer and plot the resulting orbits. 



In Figure [L3] some orbits are plotted. As expected, for fx = 0.9 one has that 
x„ — )■ X* = 0, whereas for /i = 2.85 the attractive fixed point is at x* ~ 0.65. It 
is interesting to notice that at /x = 3.2 there are no fixed points, but x„ oscillates 
between two points. We will say in this case that there is an attractive 2-cycle. 
When /i is further increased to /x = 3.5, the motion oscillates between four points- 
an attractive 4-cycle. Finally, for very large values of such as /x = 3.9, there is no 
apparent pattern for x„; indeed it will turn out that there the motion is chaotic. 
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Figure 1.4: A period-doubling bifurcation. 



It is worth investigating in more detail what happens when passing = 3. 
First, let us remark that when has a two-cycle, {f^Y = ffi ° ffj. niust have two 
fixed points: 

(1.37) f^{xi) = X2 , fi,{,X2)=Xi (f^)'^(xi) = Xi, {ft,f{x2)=X2. 

Let us now set = 3 — e. Then x* = is a stable fixed point with slope 
= s — 1. Correspondingly, the composition o has also a fixed point 
there, with slope (/^ o f^)\x*) = [e — 1)^ < 1, again stable. 

Let now /x = ?)+e. Then x* = is unstable, with slope /^(x*) = —e — 1 < —1. The 
same fixed point for the composition has then slope {f^of^y(^x*) = {e+lY > 1, and 
it is also unstable. Furthermore, as depicted in Figure [T]4| by continuity a couple of 
fixed points are created to the left and to the right of x*, and it is easy to see that 
they are stable. 

The above reasoning seems to be applicable not only when going from period- 
one to period-two, but every time we double the period of the attractive cycle. We 
have already seen that for larger /i there exists an attractive four-cycle. It is worth 
plotting the bifurcation diagram for the logistic map, that indicates for any /i the 



set to which the motion is attracted. Looking at Figure L5 we se that at several 
points fio, fii, . . . , fin, ■ ■ ■ a bifurcation occurs, where the numebr of attractive points 
doubles (i.e. one goes from a 2"-cycle to a 2"'"'"^ one). As we discussed, /iq = 3, and 
one can see from the plot that fii ^ 3.45, fi2 ~ 3.55, etc. 

Furthermore, these points seem to accumulate to some fi^o ~ 3.6 after which 
the trajectory jumps wildly between many points. To better understand this, let 



us look at Figure 1.6 One sees that the Lyapunov exponent 7(yu) is smaller than 
zero in presence of a 2"-cycle, and vanishing at the bifurcation point as it should. 
However, for /i > fi^c one sees that 7(/i) > which signals the beginning of chaotic 
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Figure 1.6: Lyapunov exponent for the logistic map for > 3. 
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Table 1.1: Values of /Xj and 6j. Table 1.2: Values of /i„, dn and a, 



behavioij^ It is also interesting to notice that between any couple of bifurcation 
points there is a point where 7(/in) = —oo, with < f^n < fJ^n- We will say that 
there the system has 2"^ -super stable cycle. 

We are not interested in what happens when /i > What is interesting, 
and we have established numerically, is that there is a sequence of period-doubling 
bifurcations at /i„ and, correspondingly, a sequence of superstable 2"-cycles at 
that converge to some jioo- 

It is not hard to compute numerically (e.g. by Newton's method) the values of 
the first few /i„'s to a good precision. These are written in Table [lT} and it is not 
hard to see that the sequence Hn — ^ /Uoo ~ 3.5699 converges, at least approximately, 
geometrically. We will call the number 

(1.38) 5 = lim ^" ~ ^""^ ^ 4.669201609 , 



10 

Clearly the same rate dictates the convergence of well. 



Feiqenhaum's 5 Notice that ^" ^" ^ is not exactly equal to 5 when n is finite. 



^For /i even larger there exists some interval where the behavior is back to periodic, with 
7(/i) < 0. These are the so-called periodic-windows, which we will not discuss here. 

^°It is interesting how Feigenbaum found this rate of convergence; he was studying the sequence 
/i„ on a pocket calculator, and needed to guess the next bifurcation point as well as he could in 
order not to waste computer time. In doing so, he realized that the convergence was geometric. 
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Figure 1.7: The dn sequence used to define a. 



Looking back at Figure 1.5 we see that the constant 6 dictates the horizontal 
scale in the sequence of bifurcations. Clearly enough, there is also a vertical scale: 
in fact, after each bifurcation, the couple of new attractive points generated that are 
generated spread out in a C-shaped figur^"] as /i increases. The size of this C shrinks 
bifurcation after bifurcation, and it makes a lot of sense to suspect that it does so 
geometrically. The "size" dn can be defined as the distance between two neighboring 
points in a superstable cycle. This can be found by looking at the intersection of 
the bifurcation sequence with the line x = Xmax = 1/2, as we will see in Exercise 



1.12 and as depicted in Figure 1.7 



In Table 1.2 the first few values of dn, as found numerically, are written. Indeed 



they converge geometrically, and one can define Feigenhaum's a as 
(1.39) a = lim ^ -2.502907875 , 

n-s>oo (i„-i-1 



which is negative because we keep into account the sign of dn, see Figure 1.7 



Exercise 1.11 (The period-doubling mechanism). Looking at Figure I.4 



we 



argued that, for the logistic map, when an equilibrium for becomes unstable, then 
two new fixed points are generated for f^. Argue that these are stable, and find the 
properties of f^ that are important for this mechanism to happen. 

Exercise 1.12 (Superstable cycles). Prove that a T^-cycle for the logistic map 

(1-40) f-^Axi)=X2, ... //i„(a;j) = Xj+i , ... /^Jxz-) = Xi , 



11 



Or U-shaped: in fact historically this goes under the name of U-sequence. 
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is superstable if and only if xj = 1/2 for some j. Argue that indeed this happens in 
between of any couple of bifurcations. 



1.2.3 Universality 

So far we have established that a certain family of maps of the interval into itself 
exhibits a sequence of period-doubling bifurcations, with geometric rate tending to 
6 ~ 4.669, leading to a chaotic behavior. What makes this story more interesting is 
that this behavior is common to many maps as well as real-world systems, i.e. it is, 
at least to some extent, universal. 

In fact, one can check numerically that the same sequence of period-doubling 
bifurcations occurs for dynamical systems defined by the maps on / = [0, 1] 

f^{x) = /i(l-x2)(2x-x^), 0</i<-|, 

16 

(1.41) f^{x) = /i sin(7rx) , < /i < 1 , 

as well as many others. What is more remarkable, the bifurcations points /i„ con- 
verge geometrically to some fi^o with a rate that tends to 5 ~ 4.669. This also 
happens also for dynamical systems of different kind, such as Mandelbrot's set, or 
Rossler's system of ODEs, see e.g. [70]. 

Furthermore, and what is more important from a physicist's perspective, the 
period-doubling cascade towards chaos occurs also in real-world systems. For in- 
stance, let us consider a fluid-dynamics experiment of Rayleigh-Benard convection, 
following Libchaber and Maurer [17] . 

Consider a box containing a fluid. The bottom of the box is kept at tempera- 
ture T, whereas the top is kept at T -|- AT. The temperature difference (or rather 
a related dimensionless quantity called Rayleigh's number) is the external param- 
eter that the experimenter may vary. When the temperature difference is small, 
heat is conducted to the colder upper surface. However, increasing the gradient, 
the familiar convective motions are generated. These consists of several counter- 



rotating cylinders that drive steadily the hotter fluid upwards,^ Further increasing 
AT leads to a more complicated dynamics of the fluid: the heat flow is not steady 
anymore, but fluctuates, as it can be seen by measuring the time-evolution of the 
local temperature at a given point in the upper surface. 

This is a discrete-time dynamical system: AT plays the role of /i, and the 
local temperature the role of $"'(x). What was found then is that, as one increases 



^^In practice, in order to obtain a stable enough convective motion, great care has to be taken in 
setting the experiment, such as picking appropriate shape and dimension of the box, and of course 
an appropriate fluid. 
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AT ~ fi, one goes from the steady temperature (fixed point) to a two-cycle, then 
to a four-cycle, and so on. Even more strikingly, the period-doubling bifurcations 
occur geometrically with rate 4.4 ±0.1 compatible with 6. This has been shown 
to occur in a number of experiments in hydrodynamics [m HHl El E3] , electronics 
E21 in ES] , charged gases flOj and chemistry [67j. 

It is worth pointing out that obtaining these experimental results is quite 
hard. On top of difficulties such as suppressing the noise and avoid generating 
chaotic behavior due to other kind of turbulence, a key obstacle is that (due to 
the geometric progression) it is in practice possible to measure only the first few 
bifurcations. On the other hand, ^""^"'^ will approach 6 only asymptotically. This 
makes these result even more remarkable. 

Of course, the word "universal" should be taken with a pinch of salt. The 



maps (1.41) are not terribly general, and indeed they share some features with the 



logistic map: 



1. ff_i{x) is regulai 
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2. ffj,{x) is unimodal, i.e. it has one maximum Xmax and satisfies 

(1.42) /^(x) > , X < Xmax & f'^,ix) < , X > Xmax ■ 

3. ffj.{x) has a quadratic maximum 

(1.43) /;(w)<o. 

It turns out that these requisites are indeed necessary, as we will also see in the 
exercises by considering some examples. More abstractly, it is reasonable to require 
some kind of regularity, since we used it in the previous section to explain period- 
doubling. As for unimodality, it is also quite clear that a very general f]i{x) with 
many maxima and minima may have a richer dynamics than the simpler examples 
we considered. As for the third requirement, it is hard to justify it a priori, and we 
will take it as an "experimental" evidence. What turns out is that maps that satisfy 
all requisites but have a maximum of higher order show the same period-doubling 
cascade, but with different universal constants in place of 6 and a, 



14 



In the next section we will see how all these features can be explained in a 
"renormalization group" framework, which also yields quantitative predictions for 6 
and a. 



^^A\\ our examples are analytic functions, which is a very strong requirement. We will not go 
into the details of how regular we need f^i{x) to be. 
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Accurate values of such constants can be found e.g. in 111) 
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1 





Figure 1.8: The tent map. 



/(l)=-a\ 



Figure 1.9: Logistic map rescaled. 



Exercise 1.13 (The tent map). Consider the "tent map" of Figure 1.8 
(1.44) 



jjx X < 1/2 

fiil-x) x>l/2 



< /i < 2. 



Prove that no period- doubling bifurcation occurs, and that in fact there is a sharp 
transition to chaotic behavior. 

Exercise 1.14 (Universahty for quartic tips). Consider the map on the inter- 
val I = [0, 1] 

(1.45) ff,{x) = fi{x - x^) {1 -2x + 2x^) , 0<fi<8, 

which, as you can verify, has a quartic maximum. Show numerically that there is 
a period- doubling cascade, and get a (rough!) estimate for the analogs of 6 and a. 
Note that this exercise can take some time. 
(Answer: 6 ~ 7.3 and a ~ —1-7) 



1.3 Renormalization group approach 

In this section we will discuss how the universal behavior of period-doubling bifur- 
cations can be explained in terms of renormalization group (RG) techniques. In 
this context, both words "renormalization" and "group" make little sense (to be 
fair, this name is a bit misleading in practically any context). They just indicate a 
rather general set of ideas of broad application, from Quantum Field Theory to the 
Physics of phase transitions. 

The reader that has some familiarity with RG in statistical Physics will find 
many similarities with what we are doing now. We will come back to this at the 
end of the section. 
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1.3.1 Heuristic 

We are concerned with several one-parameter families of maps of the interval 
/ into itself that exhibit similar period-doubling cascade of bifurcations. It makes 
sense to consider the space U of all maps with "good properties" , 

(1.46) U = {f : I ^ I , f regular, unimodal and with quadratic tip} . 

Observe that the families {/^} C U are curves in U. Clearly W is a subset of some 
space of functions, and to proceed rigorously further formalization (e.g. on the 
metric of this space) would be needed; here we will be qualitative. 

Let us chose explicitly /. Earlier we picked / = [0, 1], but to better keep track 
of the maximum of f^{x) we will set / = [— 1, 1] in such a way that the maximum 
is in a; = and takes value /^(O) = 1; we further restrict to even maps to simplify 



figures and discussion. A typical map is shown in Figure 1.9 

To single out the fundamental characteristics of the universal behavior, it is 
easier to think in terms of superstable maps: we have that, for any n = 1,2, . . . at 
/i = /i„ there exists a superstable map of period 2", and of characteristic size dn- 
The convergence of /i„ — )■ /ioo has universal rate S, whereas the one of (i„ — )■ has 
rate a. 

Therefore, universality should emerge as a property of U under the action of 
some "renormalization" operator 91 :U ^ U. This operator 

• Should relate maps with a 2^-cycle to maps with a 2^~^ cycle, 

• Should relate superstable maps to superstable maps, up to a rescaling of a. 

The first property suggests that it must be 9^(/) ^ f o f, which however is not an 
operator on U. In fact, it is easy to see that if f{x) is unimodal, / o /(x) is not. It 
is clear that some kind of rescaling is needed. 

A generic map / G W has its minimum on / exactly at the boundaries of the 



interval. Looking at Figure |1.9| define 

(1.47) a = -/(l), b = f{a). 

With a drawing, one can convince himself that the following inclusions hold 

(1.48) /([-1, 1]) = [-a, 1] , /([-a, a]) = [b, 1] , /([6, 1]) = [-a, /(6)] C [-a, a] 
provided that it is 

(1.49) 0<a<6<l, f{b)<a. 
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Figure 1.10: The renormalization operator DK. 



In this case therefore one has that 

(1.50) fof: [-a, a] [-a, a] . 

Then, there is no problem to act with on /, if we rescale all variables in an 
appropriate way (see Figure 1.10): 



1.51) 



^H(/) (x) = — /o/(-ax), 
a 



provided that (1.49) holds. In other words, (1.49) are conditions of / G W to be in 
the domain of 9i or, as sometimes it is said, for / to be renormalizable. 



Even if in general it is not immediate to see whether also 1H(/) is in the 
domain of for superstable maps things are easier. In fact, superstable maps are 
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renormalizable, and if / is a period-2" superstable map, then is a period-2" ^ 



superstable map, see Exercise |1.12[ The values /i = /ioo, where each family of maps 
gets to the onset of chaos, identify maps that have 2°°-period, and are therefore 
infinitely renormalizable. 

It is now the time to change point of view, and start studying a new dynamical 
system, where M. = U and $ = here we want to understand the orbits of the 
"points" f eU (that are, in fact, functions in an infinite-dimensional space) under 
action of This analysis may be more mathematically complicated and subtle 
than the previous ones, but we will reason by analogy with what we have seen until 
now. 

Unfortunately, it is very hard to establish the properties of this infinite di- 
mensional dynamical system, but a lot of progress can be made if we accept some 
conjectures, originally put forward by Feigenbaum [271 [2H1 129], see also [301 122] . 



1. There is a fixed point 0* G U, i.e. 9^(0*) = 0*. 

2. The fixed point is hyperbolic, meaning that the derivative of the renormaliza- 
tion operator R = c/fH at 0* has no eigenvalue of modulus one. 

3. Only one of its eigenvalues has modulus larger than one; we will call it 6. 



Under these assumptions, it is reasonable to assume that there exists an un- 
stable manifold W" of dimension one which generalizes the eigenspace relative to 
S, and a stable manifold of codimension one, that generalizes the eigenspace of 
stable eigenvectors. Let us make an additional assumption. 

4. Let S„ be the manifolds of period-2" superstable maps in U. Then inter- 
sects El transversally at 0q G W" fl Si. 

We have already remarked that sends superstable maps into superstable 
maps. Therefore we have the inclusions 

(1.52) ^"(S„+i) C Si , 

and it is not hard to imagine that all the S„'s will intersect transversally W", at 
points 4>n-i- It is also clear that the sequence 0* converges geometrically to the 
fixed point 0*, with rate 6. One can imagine that not only the points 0* — )■ 0*, 
but also the manifolds S„ accumulate toward the unstable manifold, and that their 
distance decreases geometrically with rate 6. 

The whole picture is summarized in Figure This also suggest how to 

explain the period-doubling cascade in a generic family of maps {f^} C U. In 
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fact, the sequence of bifurcations occurring at /i^, or equivalently the sequence of 
superstable maps at fin can be described in terms of the behavior of the sequence 
of manifolds S„. We will return on this later, in order to make the relation with 6 
more quantitative. 

As remarked, this whole discussion has been very qualitative. A rigorous treat- 
ment would bring us too far from the points of our discussions; the interested reader 
is invited to consult e.g. [18]. Here it is worth mentioning that, once Feigenbaum's 
conjectures are accepted, it is not hard to prove that the scenario we described 
happens. What is much harder is to establish the existence of the hyperbolic fixed 
point. Remarkably, all this could be done rigorously [171 [191 HU [23l [I5] . Finally, let 
us stress that had we relaxed the condition that our maps are even, we would still 
have found a single hyperbolic fixed point 0*, which turns out to be even. 

Exercise 1.15 (Renormalization operator and logistic map). Consider the fam- 
ily of maps on I = [—1,1] 

(1.53) f^{x) = l-fix\ 0</i<2, 

which is a suitable rescaling of the logistic map. Work out what the renormalizability 



conditions (1.49) imply in terms of jj,, and comment on what it means. 



Exercise 1.16 (Domain of the renormalization operator). We implicitly as- 
sumed that superstable maps are in the domain of the renormalization operator D^, 
.e. satisfy ( 1.4^ - Prove it. 



I. 



Exercise 1.17 (Linearizing the renormalization operator). Compute formally 
the (Frechet) derivative R = Dy{ of 91 at (p, recalling that this can be found from 
expanding formally to first order f = (p + ap, showing that it is 

(1.54) R<^ ■ (p{x) = — </?((/)(— ax)) {p{—ax) 0'(0(— ax)) . 

ct ct 
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1.3.2 Predictions from renormalization group 

Now that we have formulated the RG ideas for maps on the interval, let us try to 
obtain some quantitative predictions out of them. 

The first step is to find some information of the RG fixed point 0*. We im- 
mediately encounter a difficulty: we must understand the role of a which appears 



in the definition (1.51). It has to do with the rescaling of the x-direction, and we 
know that for the fixed point this amounts to shrinking by (negative) a. This leads 
to the identification, at the fixed point 

(1.55) a = -l/a. 

This still leaves a undetermined. However, it is easy to see how this is fixed 
by the normalization of the maximum to 0*(O) = 1. Let us consider the ansatz for 
a symmetric 0* (x) 

N 



:i.56) 0*(a;) = 1 + J]c„x2" + 0(a;2^+2)_ 



n=l 



Plugging this into the fixed point equation 

(1.57) 0*(x) = ^H(0*)(x) = a0*(0*(x/a)) , 
we find e.g. for = 3 the solution 

(1.58) a ^-2.479, ci^ -1.522, cs ^ 0.073 , cg « 0.046 . 

Going to = 6 yields a ~ —2.502897, an estimate which turns out to be correct 
up to order ~ 10~^. Incidentally, the above procedure illustrates the importance of 
the order of the maximum, which is a crucial ingredient in our ansatz: a different 
choice would have lead to a different result for 0*(x). 

Let us now try to find a more quantitative relation between 5 and the sequence 
/^n /^oo for a family of maps {/^} C W. It is convenient to denote introduce the 
short-hands 

(1.59) F(x)^/^^(x), ^(x)-^^' 



Then a function can be written, when /i is close to /ioo, as 
(1.60) f^{x) ^ F{x) + {^^- /ioo) <^(x) , 
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and similarly we can expand 

(1.61) ~ HF)ix) + (/i - /ioo) Rf ■ ^{x) . 

Recall that F lies exactly in the intersection fl {ffi}- Therefore, ^{F) G 

is closer to (p* than F, and indeed due to the geometric convergence on the stable 

manifold, ^R"^{F) ^ (p* after a few iterations. We can thus write 

(1.62) ^"(/M)(a^) ^ + (/i - /ioo) R^. ■ Vix) . 
Let us expand on a basis of eigenfunctions of R^* , 

(1.63) ip{x) = csipsix) + ^Cj(fj{x) , 

j 

where we have distinguished the eigenfunction pertaining to 6. Since all eigenvalues 
except 6 have modulus smaller than one, their eigenvectors ipj, j ^ S are sent to 
zero by R^. . Only the eigenvector ips of 6 plays a role, so that we can eventually 
write 

(1.64) ^"(/M)(a;) ~ <f>*i^) + (/i - /^oo) cs 5" <^5(a;) . 

Let us now specialize the above expression to the case fj, = fin, i.e. the case 
where the map is superstable of period 2"', and evaluate it at x = 0. On the one 
hand, we have 

(1-65) ^"(/aJ(0) = (/aJ'"(0) = 0, 

due to the presence of the 2"-cycle and the fact that a; = is a point of the cycle. 
On the other hand we have 

(1.66) (f)*{0) + {fin- fioo)cs6'^l^s{0) = - + {ftn- f^oo)S'^ Cs Kg , 

a 

where we emphasized that ks = ^ps{0) does not depend on n. Therefore, at least up 
to higher order terms in /i„ — /Xoo it must be 

(1.67) (/in -/ioo) 5" ~ — = const. , 

a Cs Ks 

for any 5, which means exactly that the rate of convergence /i„ — )■ oo is 5. 

The only thing that remains to do is to compute 5, using the expression for the 



differential R found in Exercise 1.17, and inserting the approximate result for (j)*[x) 
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found by means of the power series expansion into (1.54). The eigenvalue equation 
reads 



(1.68) 



aip5[(j){x/a)) + a(ps{x/a) (j)' [(j){x / a)) = 6 (ps{x) 



and can be solved approximately by using an ansatz for fs{x) too. The result, 
taking = 6 in the ansatze, is 

(1.69) 5^ 4.66914, 

with an accuracy of order ~ 10^^ with respect to the known result [TT] . 



We conclude this first part with some comments on what we have seen. We 
have considered systems that are described by iterated maps on the interval (that 
is among the simplest dynamics one can imagine) where the experimenter is able 
to tune one parameter fi; some class of these systems exhibit similar properties as 
the parameter approaches a critical value /ioo- We have explained this by using the 
properties of the renormalization of operator 9^. 

What is the physical interpretation of !H? When we are looking at the dynam- 
ics, yi operates a rescaling of the time-scale (/ / o /) together with a rescahng of 
the x-scale. In this sense it is similar to Kadanoff 's coarse-graining transformation 
1^ 112] : acting with the renormalization operators corresponds to changing the de- 
scription of the problem, "zooming out" in such a way as to preserve the interesting 
(universal) properties of the dynamics. 

In the language of statistical physics, we would say that /i plays the role of 
some adjustable "knob" (temperature, external magnetic field, etc.), and that at 
/ioo a phase transition occurs. The divergence of the correlation length in our case 
is mimicked by the appearance of an infinite-period cycle. All the systems at the 
phase transitions are points on the stable manifold W*, and due to this they are 
very similar. In fact, under the action of 91, all these points get to the fixed point (p*, 
so that the properties invariant under 5^ (the large scale properties, in a statistical 
system) are common to all of them. We did not investigate at all what happens to 
our maps at fioo (to avoid the complications of chaotic systems), but it is indeed 
possible to single out several universal properties. 

The role played by 6 is that of a relevant eigenvalue, and (ps{x) is a relevant 
direction at the fixed point. This means that if we perturb 0*(x) by something 
proportional to (ps{x), acting with 9^ will drive us away from the fixed point. No- 
tice that in Figure |1.11| the stable manifold separates U in two regions. Maps on 
either side of the W will have different behavior under 91 and therefore different 
properties. We have seen that all the maps "below" the stable manifold (/x < /ioo) 
are periodic, with negative Lyapunov exponent, whereas indeed maps "above" 
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will yield chaotic behavior. Again this is in analogy with what happens in statistical 
systems, where the position in parameter space of a theory determines its large-scale 
properties. 

Exercise 1.18 (Renormalization group for quartic tips). Work out, on a com- 
puter, the computation of a and 5 for maps with quartic tip. 
(Answer: a ^ -1.69030297 and S ^ 7.28486622) 



Part 2 



Functional renormalization group 
equations 

In this part we will consider a renormalization group technique for Quantum Field 
Theory (QFT) that is referred to as Functional Renormalization Group Equation 
(FRGE) or, quite improperly, as Exact Renormalization Group Equations (ERGE). 
Before doing so, we will recap some textbook notions about QFT, see e.g. [SH 175] . 

2.1 Elements of quantum field theory 

Quantum Field Theory (QFT) arose as an attempt to describe the quantum theory of 
elementary point-like particles so as to be consistent with the requirements of Special 
Relativity. One consider the Hilbert space T-L of physical states, which bears a unitary 
representation of the Poincare group, under which only the vacuum state \vo) is left 
invariant and such that its generators have positive norm. In this language fields 
are distributions which take values in some space of operators on %. Furthermore, 
since no signal can propagate faster than light, we require that fields smeared by 
test functions with disjoint supports whose distance is space-like must commute, or 
anti-commute according to the statistics they obey. 

Let us consider a free scalar field in four-dimensional Minkowski space. The 

state space is a Fock space 

oo 

(2.1) 'H = @H^, 

n=0 

where Hn is the n-particle Hilbert space (the symmetrized n-fold tensor product 
of the single particle one), and Hq is generated by the vacuum \vo). In this case. 
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Hi is generated by the complete orthonormal set of positive energy solutions fi of 
Klein-Gordon equatioiiQ 

(2.2) [r]^'^d,d, + m^) /(x) = 0, 
and the scalar product is 

(2.3) (/., /,) = I dM f:ip) f,ip) = 5.,, 
where 

(2.4) fix) = I dM fip) e-'^-\ df,{p) = ^27r5(p2 - my{po). 

The free scalar field 4>{x) also satisfies Klein-Gordon equation, and can be 
written in terms of creation and annihilation operators a\ a, 

(2.5) <P{x) = [ dfi{p) {a{p)e-'P-'^ + a+(p)e*^-") , 



which satisfy canonical commutation relations: when a^, a are smeared by a solution 
/ of Klein-Gordon's equation, i.e. a,/ = J dfifa, one has 



(2.6) F/.' 4J = = 0' 4J = 

The physical content of the theory, in this simple case, is given by the vacuum 



expectation value of any (suitably ordered) product of field operators (2.5), which 



can be explicitly expressed in terms of Green's functions, using (2.6) and Wick's 
Theorem. 

The theory of free fields is of little interest for Physics. We want to consider 
theories defined by Lagrangian densities of the form e.g. 

(2.7) L = Lfree + Lint = \ {v^'O^.^d,^ - m^') + Lint, 

where the free Lagrangian L^ee alone would give rise to ( |2.2[ ) as its variational equa- 
tion^ To describe the Physics of interacting QFTs, one considers two asymptotic 
state spaces "Hm and Tiout, both isomorphic to the space of free particles, related 
by a unitary operator, the "S'-matrix", which contains the information about the 
scattering. In this language, the S'-matrix for a free theory will just be the identity. 



^We take the metric to have signature (H ). 

^Notice that one would find modified Klein-Gordon equations for the interacting theory. How- 
ever those would involve products of 0, that is, ill-defined products of distributions, which is a first 
hint of the presence of ultraviolet pathologies. 
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The predictions of a QFT, such as scattering cross sections, will depend on its 
S'-matrix. This, under suitable asymptotic conditions, can be expressed in terms of 
the time-ordered ra-points correlation functions of the theory, by virtue of Lehmann- 
Symanzik-Zimmermann (LSZ) formula. Such quantities, in the path integral formal- 
ism invented by Feynman, are the momenta of an (ill-defined) functional measure 
over the "paths", or rather histories, of the fields configurations: 

(2.8) (0(a;i) ■ ■ ■ <P{x,)) = ^ J 1^<P {H^i) " " " 4>{xk) e^/'^'^^[<^(^)l) , 

where M = JVcpe'-l'^ is an (infinite) normalization constant, such that (1) = 1. 



The issue in (2.8) is that field configurations are weighted by a complex exponential; 
furthermore the action in Minkowski space involves hyperbolic differential operators, 
which are rather unpleasant to deal with. A much better setting would be provided 



by the "Wick rotation" (the analytic continuation) of (2.8), which formally formally 



yields the "Euclidean" theory, expressed by the Feynman-Kac formula 

(2.9) (0(a;i)---0(xfc)) = -^ j P0(0(xi)---0(a;,)e-^'^'^^[<^(^)]) . 

The Euclidean correlation functions are called Schwinger functions and, remarkably, 
they can be rigorously related to the correlations in the physical, Minkowskian 
theory [57]. In principle one can construct the QFT in a Euclidean formalism and 
reconstruct it in Minkowski space. For the analytic continuation to be possible and 
lead to a physically meaningful relativistic theory, the Schwinger functions have to 
satisfy a number of conditions, or axioms [7T| [82] . However, even Euclidean QFTs 
are far from being well defined. A number of pathologies arise when one tries to 



evaluate (2.9) to extract physical predictions. In what follows, we will always deal 



with Euclidean theories and forget about the reconstruction problem. 



2.1.1 Generating functionals 

Let us consider the Euclidean theory of a free scalar field; this amounts to studying 
the (normalized) functional measure formally defined as 

(2.10) dfici<f>) = rii ^ p^g-|/d*x<A(x)(-AW)0(x) 

where we have written to remember that the measure depends on a certain 
operator C = [—A + m^]~^j^ Since the integral appearing as a normalization is 



■^A motivation to introduce dfiQ is that one can rigorously define it as the conditional Wiener 
measure of covariance C on the space of "paths" 0, whereas Vcj) is just a formal measure. 
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Gaussian and can be performed exactly, we find that the normahzation constant is 
given by / P0e~2 /''^(~'^+"' = (det[— A + m^])"^/^. We want now to compute the 
quantities 

(2.11) {(Pixi) ■ ■ ■ (pixk)) = j due (0(a;i) ■ ■ ■ (pixk)) , 
and this can be done by introducing the generating functional Z[J] 

(2.12) Z[J] = y"rf/ice^'^'^-^(^)*(^\ 



such that 
(2.13) 



{(p{xi) ■ ■ ■ (p{Xk)) 



5 



SJ{xk) 6J{xi) 



log Z[J] 



j=o 



As is well known, the generating functional can be evaluated exactly in this 
case, by performing the Gaussian integral (by a simple change of variables), yielding 



(2.14) 



Z[J] 



62 



yd^xd%J(x)Cix-y)Jiy) 



C{z) 



p2 _|_ 



from which it follows that the two-points Schwinger function is just the free Eu- 
clidean propagatoir} Similarly, we can find the explicit expression for all correlation 



(2.11), which can be expressed by the familiar Feynman diagrams. We could intro- 



duce a different functional of the sources Vr[J], defined by 



(2.15) 



-W[J] 



Z[J], 



so that we do not have to take the logarithm in (2.13). This can be proven to 



generate all the connected correlation functions, see e.g. [90j §6. For the free theory, 
the only non-vanishing one is the two-points function. Let us introduce another 
generating functional, which will be useful later on. Observe that in the free case 



we have, using (2.14), 
SZ[J] 



(2.16) 



SJ{x) 



d^yC{x~y)J{y) Z[J], 



SW[J] 
6J{x) 



d^yC{x - y)J{y), 



so that if we define the effective field as the expectation value ^{x) = {(p{x)) 

(-A + m'^)ip{x) = J{x). 



^^^■^^ I J, this satisfies the equations of motion 



SJ{x) 

(2.17) 



^Of course, if we wanted to consider the Minkowskian theory, we would have to worry about a 
prescription to obtain the correct Green's function by analytic continuation. 
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We can invert the above relation to obtain J[v5], and substitute it in the generating 
functional This amounts to defining the new functional T[ip] = W[J[ip]], 

which, in a generic interacting theory, can be done by Legendre transform 

(2.18) T[^] = iiif (w[J] + j d^xJ{x)^{x)^ , 

from which we find the relation 

This is a constraint for the effective field ^p{x) which in absence of currents must 
solve a variational equation similar to the one for S[(j)\ in the classical theory. For 
this reason, we can consider T[ip\ as an quantum effective action. We have two 
complementary interpretations for F: on the one hand, in a diagrammatic approach, 
it can be proven that F can be obtained from the one particle irreducibl^ (IPI) 
Feynman diagrams: 

(2.20) F[(^] = I Vcj) e-^^'^+^\ 

Jipi 

where by this formal notation (borrowed from [79]) we indicate that the path integral 
is restricted in such a way that, when expanded in a perturbative series, only IPI 
graphs appear]^ We can also think of F as given by an infinite series in whose 
coefficients depend on the loop integral^ Such an expansion in effective vertices 
will be "semi- local" , since we are not just summing tree diagrams: 

(2.21) = 5^-^ /" dV---rfVr(")[xi,...,a;^]<^(xi)---<^(x„). 

n ' 

This gives a second interpretation: the full quantum theory generated by the action 
5" is equivalent to the classical (tree-level) theory for an action taken to be Smw = T, 
computed in terms of the effective vertices. As a result we can write 

(2.22) T[ip] = S[{p] + quantum corrections. 
For more on this, see §16 in [79] and §8 in [75] . 



^Recall that a Feynman graph is one particle irreducible if it cannot be disconnected by cutting 
any single internal line. 

^This was actually the way the effective action was originally defined ,36^ , whereas the functional 



definition (2.181 was given in |25l[30]; see also [39] for a derivation of the n-loop expansion of F. 
^Actually, on the renormalized loop integrals, as it will be clearer later. 
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If we go back to the free theory, we can exphcitly substitute (2.17) in the 
definition of W[J] to find that 



(2.23) 



d^x-ip{x) (—A + m^) ip{x) 



This once again shows that a free theory, even if we take into account all possible 
quantum corrections, just describes the propagation of noninteracting particles, as 
was already clear from (2.14), and the expansion (2.21) is just given by the two 

6^W[J] 



point function. Furthermore, we see that 



(2.24) 



[5J{x) SJ{y) 



which, as can be checked using (2.19), is true also in a generic theory. Furthermore, 
one also obtains from (2.18) tha^^] 

^ J Vp exp!^-S[p + yD] + j d^x^^^ 



p{x) 



yd'^xd%^{x)C(x-y)^{y) / ^^j g-5,„t [p]+/d*xp(x)(- A+m2)-p(x) ^ 



exp(-r[v9]) 
(2.25) 

that can be taken as an alternative definition of T[ip], where we replace the current 
J{x) by (—A + m'^)ip{x). The insertion of the inverse propagator has the effect of 
"amputating" the external propagators, see §2.5 in [64J. 



Exercise 2.1 (The F and W functionals). Check (2.24) for an interacting the- 
ory. 

Exercise 2.2 (The saddle point approximation). Consider the interacting the- 
ory 

(2.26) S[^]= [d'x ^ , ... , 



and work out the saddle-point approximation ofT[ip] using (2.25). 



2.1.2 Perturbative expansion of 



As we know, the difficulties and the wonders lie in the interacting theories. We will 
consider a specific theory for our examples, which is the drosophila of QFT: the A^"^ 
theory. Its Euclidean action is 

A , 



(2.27) 



d X 



:-A + m^)</) + 



4! 



^For a proof of this last expression, see [73] §18, or [BJ §2.5 where a different strategy is used: 
the above formula is taken as definition, and the familiar properties of F are derived from it. 
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where A > 0. Its generating functional will have the form Z = j rf/ic-e""^'"*. For such 
an interacting theory, we will not be able to exactly evaluate the path integral in 
general. Instead, we can perform a perturbative expansion, i.e. a formal expansion 
in powers of A. Supposing that A ^ 1, we have 



(2.28) = dfic 



^ ^ d^xct>{xY + ( I d^xct>{xY^ 



A\ J ' ' (4!)22! 

In the same way, Z[J] is expanded as a formal series in A: 
(2.29) Z[J] = e"^^ [ dfice^-^''' = g-^^e 



and the Schwinger functions are computed order by order in the coupling parameter. 
However several difficulties arise: 

1. This procedure can work only when A is small, which makes it very hard to 
address strong-coupling issues. 

2. It is not clear how small the coupling has to be. In fact it is not even clear if 
there is any nonzero radius of convergence (even in some weak sense such as 
Borel summation) for the perturbative series. 

3. As we will see, not even the individual terms in the series are well defined, as 
many of them are given by divergent integrals. 

We will focus only on the last point. The techniques we will discuss are suitable for 
perturbative renormalization, even to all orders; if we wanted to address the issue 
of convergence of the perturbative series we would need more refined tools. 

It is simple (see e.g. [59J) to derive the first terms of the expansion of, say, the 
two-points Schwinger functioi]|^ in momentum space: 



d^q A/2 

p2 _|_ jj^2 ^p2 _|_ jYi^y J (27r)^ q'^ + 711? 



(2.30) / {mm) e-'^- = - ^ / I . . + 0{\'). 



The latter term corresponds to the well known tadpole graph (Figure 2.1(b)), that 
diverges for large internal momenta q. Let 

(2.31) T ' "'"^ 



2 (27r)4 Jq g2 + ^2 



^Should we consider the amputated two-points function, we would not find the factor of {p^ 
TO^)~^ in front of the integral. 
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Figure 2.1: One particle irreducible two-points graphs for Xcp"^, up to two loops. 



This divergent expression can be regularized for instance by means of dimensional 
regularization. The result is then a function of the dimension D 

2-D/2 

(2.32) T{D) = i ^] r ( 1 



2(47r) 



D 

~2 



where m, A have been rescaled so that they are dimensionless in D dimensions, and 
yU is a mass scale. By treating the dimension as a complex variable, it is easy to 
isolate a divergent contribution as D — )■ 4, which is just given by a simple pole. A 



similar calculation can be done for the bubble (Figure 2.2), which instead gives 

A2 



(2.33) B{D) = 11 



A-D 



327r2 



E 



A-D 



Ait 

2 + ^(1) + In — 




Am 

T 



In 



1 + ^ + 1 



-| _|_ 4Th^ 1 



+ OiD-A) 



The variables s, t, u are Mandelstam's invariant momenta, scaled to be dimension- 
less: s = s/ fi"^, etc. The finite part of the diagram has also a complicated dependence 
on the momenta entering the loop (which did not happen for the tadpole due to 
momentum conservation). 



As it is well know, the trick now is that the original action (2.27) can be 



modified by adding term that cancel the divergent contributions of (2.32) and (2.33) 
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Figure 2.2: The bubble is a one-loop divergent graph with four external legs. 



at D — )■ 4. This simply amounts to adding counterterms that are the quadratic and 
and quartic in and therefore just add up to the already existing couplings. 

Now, given an arbitrary divergent graph, we can consider its one particle 
irreducible (IPI) components, in which the momentum integration factorizes, and 
renormalize each of the divergent ones. Of course the above trick could not work if 
we were dealing with, for instance, a six-legs IPI divergent graph, since there is no 
six-points term in the original Lagrangian. Fortunately, as is well known, this never 
happens in the four- dimensional theory, as can be established by counting the 
superficial degree of divergence. 

However, it is still not clear what happens when we consider a generic IPI 
diagram containing several divergent sub-diagrams. Is it still possible to factor the 
divergent contributions and absorb them in a redefinition of the couplings, or are 
new terms generated? Already at two loops there is a divergent IPI graph whose 



internal bubbles share a line (Figure 2.1(c)); allowing more loops, the situation gets 
worse. Fortunately, it is possible to show that such overlapping divergences cancel 
out at any order, and that the theory can be indeed regularized by modifying only 
the couplings appearing in the original lagrangian. The proof of this remarkable fact 
relies on ideas due to Bogolubov, Parasiuk, Hepp and Zimmermann [H |37l EHl |89] , 
and is beyond the scope of this discussion. 



Exercise 2.3 (Dimensional regularization) . 

ularization. 



Compute (2.32) in dimensional reg- 



2.1.3 Renormalization and anomalous scaling 

The previous discussion hinted that in four dimensions is perturbatively renor- 
malizable to all orders. With this in mind, let us go back to the idea of redefining 
the couplings. We use as an input in the path integral the bare Lagrangian L^are- 
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We have 

111 

(2.34) Lbare = -Zbare(t>{-^)(t> + ^T^lare^bare(t>^ + |j Afeare -^^are^^ > 

where we have introduced a wave-function renormahzation such that (t)bare — 

1 /2 

^^hare- "^^^ Lagraugiau is not physical, i.e. its couplings cannot be determined 
by any experiment. However, we want it to satisfy 

(2.35) L,,re = ^ (1 + P)<t>{-^)(t> + (1 + Q)0' + (1 + R)^' ' 

so that the physical quantities, the n-points Schwinger functions, are finite. Here 
P, (5, R are functions— counterterms— of the dimension displaying the correct pole 
structure as D — > 4. While this makes the observables finite, there is an ambiguity 
in fixing the exact value of the physical couplings, which comes from putting a finite, 
regular (in D) part in the counterterms. This has to be done through measuring 
them in an experiment, but we notice form the explicit expressions of T{D), B[D) 
that the prediction will depend (logarithmically) on the scale /i: we have to measure 
the couplings at a given energy scale and set conveniently the bare quantities. Then 
we have to expect the physical couplings to depend on the scale. 

Let us consider the n-points IPI correlation function F*^"); observing that 

1/2 

(f^bare — -^bare^i ^^^^ rewriting Z^^are function of as we have 

(2.36) r(:i,(pi, . . . ,p„;mLe, >^bare) = ^'"/^r^") (pi, . . . , p.; (/.) , A(/.)) . 

Notice that the right hand side depends only on physical coupling, and since the 
l.h.s. is independent of /x, ^^^f^ — 0, we have the equation for the physical quantities 

f d dm? d dX d ndlnZ\ ^f„^,^ r,, , 

The same is true if we consider logarithmic derivatives. Call e = 4 — D, and define 

dX 

(2.38) l3x{X,m^,e)=fi—, 

(2.39) ^^(A,m^£)=/.^^, 

(2.40) 77^(A,m^£)=//-^, 

recalling that A = A only when e — )■ 0. Remark that once the l.h.s. (the (^-functions) 
are known, these equations determine the evolution of the couplings. 
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Let us introduce a dummy parameter s to count the powers of momenta appearing 
in r'^"^ p^sp. On dimensional grounds 

(2.41) + 2m2^ + s£ - [A-n+e{n - 2)]^ T^-\sp; /i, m^fi), A(/i)) = . 

Putting everything together, we have the scahng equation (in the hmit e — > 0) 

+4-n^ r(")(sp;/i,m^A) =0, 



supplemented by (2.38 2.40). This is Callan-Symanzik equation. The interpretation 



is that the same n-points function, considered at different energy scales, changes 
only by a rescaling. In a classical theory, we could have written immediately the 
solution: 

(2.43) r(") {sp, m. A, /i) = s^'™^ T^") {p, A, m, /i) , 

Remarkably a similar scaling holds in the quantum case. Once we find the /3- 



functions perturbatively from the counterterms by imposing the equality of (2.34) 



and (2.35) [71l[76|, we can integrate (2.38), and finally use them to solve (2.42) by 
the method of characteristics (e.g. [20j| §7.3), yielding the scaling 



(2.44) 



^(")(sp,m,A,/i) = s^-"e-t/i'^*^^^ r(")(p, A(s), m(s), /x) 



where the extra exponent gives the anomalous dimension, as is immediate to see 
in the simple case when r/^ is constant. Observe that, for a scalar theory in four 
dimension, the presence of an anomalous scaling is particularly relevant for the 
4-points function, which has zero classical dimension. 



2.1.4 Renormalizable and non-renormalizable theories 

Let us briefly recall a criterion for perturbative renormalizability, which we men- 
tioned to justify the absence of six- legged IPI divergent graphs in A0^. Consider a 
scalar theory in arbitrary dimension D, of the form 

and a IPI graph with L loops, V vertices, / internal lines and E external lines. Let 
the superficial degree of divergence 5 of a graph be given by D (from the integral 



(2.45) 
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/ d^p) times the number of loops, minus two (in the scalar case) times the number 
of propagators appearing in the integration. It is not hard to see that then 

(2.46) 5 = D - ^{D -2)E + V (^^^^^^D - N 

It can happen that a graph converges better than what its superficial degree 



of divergence would suggest^ However, by a theorem of Weinberg, ii S < the 
graph is convergent. Therefore we can give the following perturbative and superficial 
classification of quantum field theories, based on proliferation of divergent graphs: 

1. Non-renormalizable, when — N > 0. 

2. Renormalizable, when ^^^D — N = 0, so that 5 does not depend on V; this 
is the familiar case of in four dimensions. 

3. Super-renormalizable if ^^^D — < 0; in this case the UV behavior gets 
better and better in higher loops diagrams, and there are only a finite number 
of divergent ones. 

4. Finite, if there are no divergent diagrams at all. 



What is special about renormalizable theories? Let us rephrase (2.46) in terms 
of the mass dimension of the couplings. The coupling g^^^ for the g^^^'cf)^ term has 
mass dimensioiJH] 

(2.47) = [^?(^)] = D- , 

and it is immediate to see that if Ajv < the theory is non-renormalizable. This 
can easily be extended to more physical theories involving fermions and gauge fields, 
see [79] §12.1. Therefore if we consider a general theory, allowing for any term 
compatible with the symmetries, we will have an action of the form 

(2.48) S = S,en + S^onren=J29R(^^[<P]+ Yl ^^nI^^^M' 

iSren jSnonren 

where = [s'nr] < 0- If '^o"^ assume that there is a common (large) energy scale 
M such that g^^^ ^ we automatically get that, when considering processes 



^''This frequently happens in presence of symmetries, an example being the light-by-light scat- 
tering in Quantum Electrodynamics. 

^"'^It is convenient to define the operators O to include the a;-integral, so that Okin ~ 
\ J d'^xd^4)d^^(t> is dimensionless. 
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of energy <^ M^, the effect of irrelevant operators on physical observables is 
suppressed by (M^/s^)^'/^ ^ 1. Therefore a first reason for the importance of 
renormalizable theories is that at sufficiently low energy they capture the leading 
physical effects to a good accuracy. Furthermore, they can be defined by measuring 
a small number of physical parameters g'^ , and hence are very predictive. Irrelevant 
operators can be taken into account as higher order effects, which however require 
fixing additional parameters [T3] . 

This does not mean that non-renormalizable theories are useless. In fact, 
when the number of renormalizable couplings is reduced by symmetries, even to 
zero, it can be important to consider the lowest order non-renormalizable ones, as it 
happens in chiral perturbation theory, which has a remarkable experimental success 
describing quark bound states (e.g. pions) at low energies ^B]. However, when we 
reach the threshold energy ~ M, any possible coupling has to be considered, and 
the theory is no longer predictive. Usually this is regarded as a breakdown of the 
theory, meaning that it has to be replaced with a more fundamental one. For the 
effective theory of pions, the fundamental theory is Quantum Chromodynamics, but 
in principle nothing forbids that the fundamental theory cannot even be expressed 
as a Quantum Field Theory. 

Exercise 2.4 (Superficial degree of divergence). Using the topological proper- 



ties of Feynman graphs for a scalar theory, prove (2.46). Drawing some examples 
will help you. 



2.1.5 Behavior at fixed points 

The two previous sections indicated that the renormalization techniques in QFT 
have similar features to the dynamical systems we described earlier. In fact, equation 



(2.48) looks a lot like the subdivision of eigenvectors of a hyperbolic fixed point into 
relevant and irrelevant, where plays the role of critical exponents, specifying how 
the corresponding operator behave as the energy scale s is varied, depending on 
their sign. 



On top of that, Callan-Symanzik equation (2.42) can be added to the previous 
qualitative discussion to specify the case where the classical dimension is zero. In 
that case it reads (if we also neglect the mass) 

where we have introduced a Z^-rescaled coupling At vanishing coupling the 
/3-function is zero. We see that the stability property of the fixed point as s in- 



12 



This actually is not important for A0^, because at small coupling rj^ is sub- leading in A. 
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A Bw 



A Bie) 




(a) f3{g) > and increasing, (b) f3{g) < and decreasing. (c) Non Gaussian FP. 



Figure 2.3: Possible qualitative behaviors of P{g). 

creases depend on the sign of the /3-function. The fixed point (FP) that we are 
talking about is obviously the free theory that as we have checked explicitly does 
not get any quantum correction, and therefore does not renormalize. Of course the 
stability analysis could in principle be done at any fixed point, but our perturbative 
techniques allow us only to deal with the free (or Gaussian) one. 

Let us comment a bit on the possible behaviors, based on the sign of the 
/3-functions. The most interesting couplings to study perturbatively are the ones 
having zero mass dimension at the Gaussian fixed point (GFP), the marginal cou- 
plings for which quantum corrections are crucial. Suppose to consider a theory with 
some renormalizable couplings and a single marginal one, g = g- To classify behav- 
ior of the theory close to the fixed point we can assume ^g = f3{g), provided that 
we restrict to energy scales much larger than the mass. We assume that the physical 
regime for the coupling is g > 0. Qualitatively, we have 



1. Figure 2.3(a) , j3{g) > 0. In this case, g{t) increases toward high energies. Even 



if we measure the coupling to be small at low energies (as happens for the fine 
structure constant acm), we have to expect a UV regime where perturbation 
theory is not be applicable. On the other hand, in the IR the theory would 
just sink in the Gaussian fixed point. Remark that if (3{g) grows with g, it is 
possible that the solution of the differential equation does not exists for any t, 
i.e. the coupling blows up to infinity in a finite time. 



Figure 2.3(b)[ (3(g) < 0. Since it is the same scenario up to one sign, the 
preceding discussion is valid up to swapping the IR and the UV limits. In this 
case, the coupling would become smaller and smaller at high energies, and the 
theory is asymptotically free, which is the case of e.g. QCD. 



3. Figure 2.3(c 



l3{g*) = 0, g* ^ 0. We have two possibilities, depending on 
the sign of /3(A). Let us take it to be positive as in figure. In this case. 
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taking an initial condition go < g* the coupling would grow at high energies, 
approaching g*. The theory would be well behaved both in the UV and the 
IR, even if, depending on how large g* is, it may not be possible to describe 
it just perturbatively. The same picture holds when f3{g) < 0, up to swapping 
IR and UV. We call these asymptotically safe theories. 

A generic dynamical system can have a much richer behavior. If we consider 
many coupled /3-functions, limit cycles or even chaotic behavior may arise. For 
certain classes of QFTs it is possible to constrain these behaviors |HZl 113] , at least 
to some extent [21] ; this is an active topic of research [SI |35l [31] , see also [5l] . 

2.2 Wetterich's FRGE for scalar fields 

Here we will introduce the Functional Renormalization Group Equation (FRGE) 
due to Wetterich [HOi |8T], in the case of a scalar field. While it is not hard to 
formally derive the FRGE, it is more involved to understand how to apply it. 

Roughly speaking, the idea is to obtain the /3-functions from the change in 
the effective action F computed up to two different scales. Rather than considering 
all quantum fluctuations, we want to integrate over the ultraviolet (UV) degrees of 
freedom down to the scale k: 



(2.50) 




which is however ill-defined in the UV, so that we must put a UV cutoff at A 
(2.51) Zk,A ~ / n ^'^(P) • 

While the above expression depends on the unphysical cutoff A, its variation with 
respect to k does not 



(2.52) 




and could be used to investigate physical properties. 

The above argument is only qualitative, since working directly on the path 
integral measure is unpractical, and we would also need to consider the presence of 
a current to generate the whole theory. However, it is possible to find a differential 
equation for the scale dependent family of functionals F^fv?], i.e. to express the 
equivalent of the right hand side of the above equation in presence of a current in 
a nice way. More extensive reviews on this subject are [TJ O [32|, |58]; interesting 
complementary reading on somewhat different approaches is [85| [H] and [61] . 
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2.2.1 Derivation of the FRGE 

In all that follows, we will deal with a single scalar field, but everything can be 
straightforwardly extended to multiplets of fields. We start from the UV cut-off 
generating functional 

(2.53) e-^-W = / PA(/)e-^[^l+-^-^ , 

A J 

where I'a^ is the functional measure with a formal UV cut-off at energy scale A, 
A/a = / r'A0e~'^f''°'=f'^^ and we introduced the short hand J ■ = / d^x J{x) (j){x). We 
obtain from it a one parameter family of generating functional 

(2.54) e-^^.^!-^ = / jy^^^-S[4>]+J-<t>-AS,[<t>] ^ 

A'fe,A J 

where now Aa;,a = / Vj>^(f)e~^''^<'^'l'^~^^''^'^^ . A5'fe[0] is a modification of the action 
which suppresses low energy modes {E < k) in the path integral. Instead of using 
a sharp cutoff in the functional measure, this can be done by setting 

(2.55) ASS] = / Hp) T^kip) H-p) = / Hp) rip'/k') H-p)- 

The regulator TZk{p) is determined by the choice of the shape function r{z). The 
requirements on r{z) are to be monotonic and to satisfy 

(2.56) r(0) > 0, lim Tlkip) > 0; 

(2.57) lim r{z) = 0, lim 7^fc(p) = 0; 

(2.58) r{z) > 0, < ^ < 1, lim Tlkip) = ^- 

fe->A->oo 

This is to ensure that TZk implements a IR cutoff, which is a sort of mass term if 
r{z) < oo; then we impose that the regulator vanishes when is in the UV; it is 
also intended that not only r{z) has to become very small when > 1, but also the 
decay to zero should be fast. Finally, when we send A; — )> A — > oo, we require the 
path integral to be dominated by the quadratic part; this can allow to set a sort of 
initial condition on the ffow equation. Particular choices of r{z) can be of the form 



(2.59) polynomial r{z) — z °', a>0, 

(2.60) exponential r(z) — — g , /3 > 1 , 

— 1 

(2.61) optimized r{z) = (1 — z)H{l — z) , 
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where H is Heaviside's step function. The last cutoff, due to Litim is called "opti- 
mized" in a sense described in [50t ISTj [52] . and we will be using it extensively. The 
different cutoffs are plotted in Figure 2A^ together with the respective shapes of the 
regularized propagators for k = 1 

1 



(2.62) 



However, provided that the regulators meet the mentioned requirements, the deriva- 
tion of Wetterich equation does not depend on details of r{z). 

From now on, we will always assume the presence of an UV cutoff A, but just 
write Wk = W^^a (and similarly A4, Tfc, etc). We consider the Legendre transform 
of Wa;[^], which is written in term of the fc-dependent mean field (p{x) = 



(2.63) 



inf 

J 



Wk[J] 



To understand what this is, we can evaluate it exactly in the free case, finding 

1 



(2.64) 
(2.65) 
(2.66) 



Wk[J] 
ffcM 



rf^xd^y J{x)Cu{x-y)J{y) 

ip-z 



(2vr) 



D 



p2 _|_ ^2 _|_ l]^2^ 



^ j d^x Lp{x) (-A + + A;V(-A/A;^)) Lp{x) 




(a) The regulator Rk{p^) for different shape (b) The (normahzed) modified propagator 
functions. Dotted hne: polynomial, a = 1. F{p'^) for different shape functions, as in the 
Light solid: exponential, /3 = 2. Dark solid: left panel. The grey thin line is the unmodified 
exponential, /3 = 1. Dashed: optimized. propagator. 



Figure 2.4: Regulator and regularized propagator for different shape functions. 
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To cancel the explicit /c-dependence from Tkl 
dependence), we define 



(which in the free case is the only 



(2.67) 



We retrieve the full quantum effective action when — )■ 0, when ASk vanishes. In 
the opposite limit, by saddle-point method we have that roo[v5] ~ S[ip]. 

With this definition, we can find the one loop effective action r^,^°°^[(y9] for a 
generic interacting theory by saddle point approximation: 



(2.68) 



tIIoop r 



lot 



(j,=tf 



which is an infrared modified version of the familiar expression and reduces to that 
when /c — )■ 0. Its differential form is inspiring: 



(2.69) 



1, 



-Tr 



k 



dUk / 6\S[4>] + ASk 
dk \ 5(f) 6(f) 



From this we can extract the one-loop /3-functions of the theory similarly to what one 
does in Wilsonian renormalization. We must expand the one-loop effective action 
on a basis of operators, assigning a coupling to each of them: 



(2.70) 



plloop 

k 



Some of them will be the ones already appearing in the action, but most of them 
will be new. Then we have to compute the trace on the r.h.s., and match the terms 
appearing in the action with the ones coming from the trace, obtaining 



(2.71) 



k 



d9k 
dk 



Pi!\9k), k 



dk 



Afgf + P'^>{gk)^P''K9k) 



~{K 13 



where in the latter expression we consider dimensionless couplings gj! 
is the mass dimension. Remarkably, this can be extended to all-loops. 



where A 

We will introduce the renormalization time t = log(A;/A), so that 4i = k 



dt 



dk- 



Differentiating (|2.54|), we have the identities 
(2.72) 



dWk[J] 



dt 



/dASk[^]\ 



J fixed 



\ dt / 



-Tr 



J 



dTZk 



j- 



dt 



^^A more physical analysis would require appropriately scaling each coupling by powers of the 
field strenght renormalization Z^/^. 
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where the trace stands for integration and, if we were to consider multiplets of 
fields, summation on internal indexes. Then, recalling that ip = (0) and using the 



straightforward generalization of (2.24) 



(2.73) 



6(p{x)6(p{y) 



5J{x)5J{y) 



where is the current used in the Legendre transform, we have 



dt 



(2.74) 



dt 
dWk[J^] 



j,JW[J]dJ^{x) 



6J{x) dt 



d^'x^ix)'^'^''^''^ 



dt 



dASkl 
dt 



dt 



dt 



dt 



-Tr 



dUi 



2 

2 
1. 



-Tr 



-Tr 



((00) J. - (0) 



dlZi 



dt 



-Tr 



dt 
6^W[J] 



SJ5J 



dlZi 



dUk 
dt 



-Tr 



dt 

-1 



dTZk 
~dt 



-Tr 



dt 



7^fc + ^(2)[ 



The above Functional Renormalization Group Equation (FRGE), due to Wetterich, 
is a differential equation for the one-parameter family of functionals F^. Its solution 
Ffc describes the flow of the effective action (hence the name flow equation) in the 
theory space under changes of the cutoff scale. The theory space is the space of all 
functionals of the fields compatible with the symmetries of the theory (the analog 
of the space U for maps on the interval). For instance, if we assume the theory of a 
single real scalar field has a symmetry for o —0 a sensible definition of a theory 
space could be 

(2.75) T = span |(— A)^0^" : p > 0,ri > 1 and permutations} . 



Some further remarks are in order: first, the above equation is formally "ex- 
act" (in the sense discussed), and in particular accounts for arbitrarily high loop 
effects. This may seem a bit surprising, because (2.74) has a one- loop structure. 
However, the equation is exact in the full theory s pace. Let us illustrate what this 
means with an example: in the theory defined by (2.75), we denote by A^"'''^''^'' the 
couplings corresponding to n powers of (j) and p laplacians (the permutations are 
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labelled by o\ Then we would have an infinite system of /3-functions 



(2.76) { |A(4'0) =/34(A(2.P'-),A(^'f''^),A(6'P.-)) 



so that even if the beta function for the y?^ coupling does not explicitly depend on 
the coupling of y?^, it does implicitly via e.g. A*^^-*. Clearly the exact solution of 



(2.76) is very difficult— actually as difficult as the exact solution of the quantum 



field theory, since we can imagine the quantum theory to be implicitly defined by 



(2.74), since F contains all the IPI diagrams. Some approximation will be needed, 
but this needs not to he perturbation theory. 

We completely lost track of the cutoff A. Indeed the flow equation does not 
contain it and appears well deflned, as the trace converges both at large and small 
momenta due to the regulators. However when we integrate it with respect to t, we 
have no guarantee that the resulting flow does not blow up. If, however, we flnd 
that the flow equation can be integrated up to arbitrary high scales, no problem 



arises. This is the case of what we have deflned in Section |2.1.5| as asymptotically 
safe theories. In that case, in T there will be two flxed points, and a certain num- 
ber of theories that flow from one to another under the RG flow. These span a 
critical manifold, whose dimension counts the physical parameters that must be 
measured to make a prediction. When the theory is deflned at the scale /cq by set- 
ting the physical values of couplings g^^\ko), . . . g^^'^ko) corresponding to operators 
Oi, . . . Omi we have to expect that at a different scale the theory is deflned in terms 
of a possibly inflnite number of operators Oj, whose couplings are determined from 

^7W(A;o),...^7("^)(fco)H 

Therefore, Wetterich equation can be useful both as an alternative method 
to flnd approximate expression for the /3-functions and as a tool to study the RG 
properties of functions that seem not to be renormalizable perturbatively. 



2.2.2 Expansions and truncations 



In practice it is not possible to address the full equation (2.74), and one has to adopt 
some approximation scheme, which generally relies on the choice of an ansatz. A 
natural way to do so is to expand it on a basis of the theory space, for instance 



taking the monomials appearing in the deflnition (2.75) as basis vector, and retain 



m^Z = 5(2,0) and Z = 



^^To make contact with the familiar notation, \Z^ 
This will become hopefully clearer in the next section. 

^^In a geometrical picture, this corresponds to saying that the integral curve € T needs not 
to be contained in any finite linear subspace of T. 
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only some of those in the ansatz, which will not be stable under the RG flow. There 
are two natural strategies: 



1. A derivative expansion. This amounts to retaining terms involving up to a 
certain number of derivatives. The rationale for doing this is that derivative 
operators make the explicit computations of /3-functions much messier, and 
often already the "local potential approximation" , which basically amounts to 
the zeroth order of the derivative expansion together with the kinetic term, 



(2.77) 



/■ r 1 °° 

d^x --Z^^Av + Y,Zl 

n=l 



X 



(2n) 



(2n)! 



2n 



can be used to make quite accurate physical predictions. Observe that we have 
replaced the coupling constants g'^^^^ — )■ Z'^X^^^\ The lowerscript k reminds 
us that the couplings are scale dependent. 

2. Expansion according to mass dimension. In this case, one fixes a mass dimen- 
sion — A < 0, and considers all the terms whose couplings have [(7] > A. For 
instance, the truncation corresponding to A = 2 in D = 4 is 



d X 



-Z^-^A^ + Z,A2'2l<^AV - Zlxf^'^^^^''^^^ 



(2.78) 



2\(4) 



"rZimi-i^' ^ ZiX^ 



4! 



+ ZiX 



3\(6) 



6! 



This can be convenient in the vicinity of the Gaussian fixed point, where the 
discarded couplings are irrelevant and should have little effect on the flow. 



In any of the above cases, however, the equation is no longer exact. In par- 
ticular, any prediction of the flow equation (fixed points, critical exponents, . . . ) 
will depend both on the truncation and on the explicit form of the regulator TZk. 
A way to make sense of the predictions is to consider larger and larger truncations, 
and different regulator schemes, and check which predictions are robust under this 
changes. Only in exceptional cases it is possible to make statements that take into 
account the flow in the whole theory space [BB], and even in such cases it is very 
hard to provide a rigorous mathematical formalization of the flow equation in T. 

Let us also mention that not all couplings are physical, since changes of variable 
in the path integral do not modify the physical prediction, but lead to redefinitions 
of the couplings. Such couplings are called redundant. 



Exercise 2.5 (FRGE for ri-points functions). We can use (2.74) to retrieve 
the behavior of the n— points functions, allowing for a graphical representation in 



50 — Part 2. Functional renormalization group equations 



Feynman diagrams. For instance, show that the two-points function is given by 



dt 



=Tr 



dUk 
~dt 



(2.79) 



-It. 

2 



k 

dTZk 
dt 



r(3) 



r(3) 



+ 



r 



(2) 



7^^ 



r{4) 



2.3 Applications of the FRGE 

We will show how the FRGE we have derived can be used to derive /3-functions, 
comparing those from the ones found by a perturbative approach, and then applying 
them to non-perturbative questions. 



2.3.1 Revisiting Xcf)^ 



Let us consider the ansatz 



(2.80) 



d^x 



1 ^ . ^ m 



,2 I -r2-^fe .A 



Then, equation (2.74) almost coincides with (2.69). The only difference is that on 



the r.h.s. of the latter we have S^'^\ whose coupling constants have no explicit 



fc-dependence. We now define the regulator in momentum space, in terms of (2.61) 
(2.81) TZkip^) = Zk - p^) H{e - p") . 

The flow equation takes the explicit form 



(2.82) 



dt 



-Tr 



2A;2 + ry(A;2 + A) 



Hi^ + A) 



where ^ = ^ log and A disappears from the denominator due to our choice of IZk ■ 
For the moment, let us not consider the running of the kinetic term, so that we can 
put r/ ^ 0. To recover the field operators of the potential, we can restrict ourselves 
to constant fields, ^p{x) = (po- In this way, the trace can be easily computed as a 
momentum integral. In fact, if Vd is a D-dimensional volume element, we have 



dt ''^ 

(2.83) 



1 
2 

V, 



d^p 
.(2^ 



2fc2 



f.D+2 



D 



(47r)^/2r(l + L'/2) k^ + mj 



k'^ + ml ' {k'^ + m\y 



1 - 



+ 
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where in the last hne we have expanded the fraction in powers of the fields]^ Notice 
how terms of the form Lpjp are generated from any n, reminding us that we are only in 
a small sector of T. From the above we find the beta functions for the dimensionless 
couplings = mlk~^, Xk = Xkk^~'^, where the explicit fc-dependence has washed 
out. 

(2.84) l'">=fe = -2'"g- (4,)or(i + |)(lAy 

which in four dimension gives 

/„ „„\ d ^ o / „ A,_9 Xh d ~ 3X^ 

(2.86) _™'-«(-2+_)™J-_, ^A.«— . 

This has to be compared with the familiar one-loop result (see e.g. [59], §4.7) from 
perturbation theory with the mass independent renormalization. There are two 
manifest differences: first, we find higher order contributions in = m^/fc^ to 
both /3-functions. Second, there is a discrepancy in the running of the mass term at 
zeroth order in m. 

The former difference can be understood as a scheme dependence due to the 
infrared cutoff, which however does not alter the qualitative behavior of the flow, 
since m^/fc^ ^ 1. In fact, it is easy to check that such depend on the choice of the 
regulator. The latter discrepancy amounts to a quadratic divergence, as it can be 
seen in terms of the dimensional quantities, k-^m\ ~ fc^/lGyr^ . Recall that these 
do not appear directly in dimensional regularization. 

We still have to extract the running of the wavefunction. For this purpose it 
is not enough to limit ourselves to constant fields. We need an x— dependent fluctu- 
ation term, such that ^p{x) = ipo + 'f>{x). Then we cannot straightforwardly perform 
the momentum integration, because there are a number of differential operators, 
acting to the right on the (p{x). However, we just need to use the commutator 

(2.87) [-A, <^(x)] = -A (^(x)) + 2id^ (^(x)) , 

and the cyclic property of the trace to sort them to one side and write formally 

(2.88) Tr \A{x) B{id)] = '^{x\A{x)\x) {x\p) {p\B{p)\p) {p\x) . 

x,p 



^^Accidentally in this case such an expansion equals one in powers of A, but the rationale for 
the expansion is to be able to match the monomials on the left and right hand side. 
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In the case of a local potential approximation for a scalar field (2.82) we can 



sort all the derivative operators to the right just by the ciclic property of the trace 
for the terms contributing to the flow of rj, i.e. for the ones quadratic the field, 
which means that it was actually consistent to set r/ = 0, consistently with what we 
know form one-loop perturbation theory. 

The higher loops effect come into play through the flow of irrelevent cou- 
plings. For instance, it is easy to see that including a term such as J —ip^Acp = 
3 / (/?^ d^ipd^^ip into the ansatz would yield a nonzero /3-function for i]. Let us sketch 
the computation. We consider 



(2.89) Tk[^] = I rf^x 



1 a(^'^) 



2 ^(4) 



To be completely consistent in D = 4, we would have to consider a more general 



action such as (2.78), but let us stick to the simplest case for simplicity. Then 



(2.90) |r.M = It. 



2e + r]{e + A) 



H{k^ + A) 



where all the derivative operators act to the right and the cyclicity of the trace is 
understood. To find the /3-functions of the potential we restrict again to constant 
configurations ipo (without setting rj = 0): 



(2.91) 



2P + r]{k'^ 



F + ml + \Z^ \t'\lp' + Ai^V^^ 



+ 0(^) 



from which, by the usual expansion, we can extract the scaling for the potential: 



(2.92) 



(4) 



(4 + D)(2 + Z}-r7) 



;(4) 



(l + m2)2 



(2.93) 



P 4 2ri)\ -4(4^)D/2r(4 + D/2) 



D(4 + L)-2r/) 



(l + m|)2 



(J>+4)(r»+6)(D+2->;) ^^(4)^2 , 

(l + m?)3 + 



^ 2£)(D+6)(D+4-77) ^(4)^(4^1) g(D+2)(g+6->7) ^ ^ (4,1) ^ ^ 
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Figure 2.5: The one-loop contribution to the wavefunction renormahzation in A0^. 



Finally we compute the two remaining /3-functions, by expanding (2.90) in 



terms of (p{x) = (po + ^{x) and using ( 2.87[ ): 

2{D + 2-7]) 



(2.94) 



(4,1) 



d r 



dt 



(4,1) 



(27r)^/2r(2 + D/2) (l + m2)2 ' 

(4,1) 



<2D-6-2r^)Ar^ + 



8(47r)^/2r(3 + D/2) 



X 



(2.95) 



X 



{1 + rhiy 
which yields a nonzero flow for t]. 



{A + D){2 + D- //_) 3^(4)3^(4,1) ^ D{A + D-v) 



{4,1)n2 
k ) 



Exercise 2.6 (Regulator dependence). Repeat the FRGE calculation to find the 
(3-functions for (2.8(^ using the more general regulator 

(2.96) RJp; a)=aZk ( k'^ - -pA H(P - -/) , a > , 

\ a J a 

to show that the quadratic divergences are non-universal. You should find 



(2.97) 
(2.98) 



d 
dt 



ml 



2ml 



327r2 [a + m 

3\2 



d~ _ 3 
dt ^ IGvr^ [a + m?Y 



3A^ 



-2 + 



A 



167r2 



Afe 



32Q;7r 



2 ' 



Exercise 2.7 (Wave-function renormalization) . Consider Xcfr' theory in D 
6. Perturbation theory calculations (JES^, §7.5^ yield 
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where the running of t] comes from the graph in Figure 2.5 . 

Using Litim's regulator, show that the P -functions and rj predicted by the FRGE are 



(2.101) 
(2.102) 



dt 



ml 



-2ml + 



V ■ 







3(47r)3 (1 + m2)3 ' 



(4vr)3(l + m2)4 



Explain the apparent discrepancy. 



2.3.2 The Wilson-Fisher fixed point in D = 3 

The techniques that we have developed up to now will allow us to consider a 
rather complicated statistical system. It is an experimental fact that several three- 
dimensional magnetic systems exhibit a second order phase transition (PT) when 
their temperature T approaches a critical value Too. In that vicinity, for a class of 
them, it is observed that the correlation length ^ diverges as 

(2.103) ^ ^ {T - ^ 

where we introduced the reduced temperature 'd = {T — Too) /Too- A typical model 
for this kind of magnets is the Ising model that, in the approximation of continuous 
spins (and in zero magnetic field), can be described by Euclidean X(f)^ theory in 



D = 3. However, the behavior (2.103) is universal, i.e. common to many magnetic 
systems, which may include more general spin interactions. 

Our experience suggests that this behavior should be explained in terms a 
property of the theory space T common to all these theories, namely the existence 
of a fixed point with one relevant eigenvalue ( (related to z/). Then there will be a 
codimension-one stable manifold, and the phase transition at T = Too will happen 
as the one-parameter curve in T described by the one-parameter family of QFTs 
under consideration intersects W"*. This will be confirmed by the RG analysis. 

In a statistical theory we are interested in the long-wavelength (IR) behavior, 
i.e. /c — )■ 0. When we integrate down from a scale ki to k2 < ki, the free energy 
will transform inhomogeneously, offsetting by a constant term which represents the 
energy of the modes that have been integrated out. Schematically, recalling that 
the free energy is given by T^^\ 

(2.104) rW(A;2; ^2) ^ AJ",,,, + (^2/^1)^°'"'' ^^^""^ r^'\k^■, , 

where the offset AJ^kiM a regular function]^ If we restrict to the singular con- 
tribution to the free energy, therefore, we get a homogeneous scaling equation. The 



17. 



This is similar to the offset ASk that we subtracted from f ^ when deriving the FRGE. 
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scaling exponent is simply given by the dimension, since F^"^ has no external lines, 



in agreement with the solution of Callan-Symanzik equation (2.44). As for the rela- 
tion between -t^i and '&2, just like what happened with iterated map of the interval, 
since we are approaching the stable manifold its scaling is given by the relevant 
eigenvalue at the fixed point, so that 

(2.105) r'il^Xs k; ^) ^ r'Sl^ ik; s< ^) , 
where s = k2/ki, or equivalently 

(2.106) T^lm/^iV^^k; ^,) ^ {A/^2f^T^Uk; ^i) . 

This indicates that, close to the PT, the relative scaling of wavenumber and tem- 
perature is A; ~ 'd~^^'^, so that the correlation length diverges as 

(2.107) ^ ( = -l/u. 

Now that we know how to relate ( to the physics of the problem, it is time 
to investigate T. As always, there is one fixed point in the theory space when the 
theory is free. Clearly it is not the one we are after (the Ising model at the PT is not 
a free theory!). There, both the 0^ and (p^ couplings are relevant, i.e. they sink into 
the GFP going to short wavelengths. This means that if we define an microscopic 
theory with given values of the corresponding couplings in the UV, these flow away 
from the perturbative region in the IR, while all n-points interactions are generated 
as effective vertices. 

To se what happens non-perturbatively we can use the FRGE [811 E]- The 



simplest way to proceed is the ansatz (2.80) which at D = 3 gives 
(2.108) 

d ^ 2 ~ 2 d ~ ~ A? 



dt " " 67r2(l + m2)2' " " n^l + miy 

Besides the Gaussian FP, there is a nontrivial one at 
(2.109) A* ^ 7.7627, m^* ^ -0.0769 . 

and the critical exponents can be estimated from the Jacobian 



(2.110) JaCNGFP 



-1.6667 -0.0198 A / -1.8426 
-25.229 1 ~ [ 1.1759 



The estimate for the irrelevant exponent gives 

(2.111) C~ -1-8426 =^ z/^ 0.5427^ 
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0,5 



N 



Figure 2.6: The critical exponent v at the Wilson-Fisher FP, in various polynomial 
truncations of degree A^. The dashed line is the prediction of seven-loops perturba- 
tion theory. 



which is not far from the measured value. 

To be sure that we are dealing with a physical effect, we need to consider more 
general truncations, and ask ourselves: 

1. Is the Wilson-Fisher fixed point present in all the truncations considered? 

2. Is it the only non gaussian fixed point? 

3. Does the dimension of the attractor at the fixed point change? 

4. Are critical exponents such as v well behaved when considering larger trunca- 
tions? 

Considering more general ansatze one sees that, even if some other fixed points may 
appear, the Wilson-Fisher one is the only one to persist in all truncations. Further- 
more, it appears that there is only ore relevant direction, and the estimate for v is 



quite stable, as shown in Figure |2^ There we compare the estimate with the result 
of resummed seven-loops perturbation theory v = 0.6304, finding a discrepancy of 
about 3% [90]. Remark that we have taken the local potential approximation where 
1] = whereas the the perturbative calculations gives rj* = 0.0335. 

We have therefore established a scenario similar to the one of Section II. 3[ 



This was first done in this context by Wilson and Fisher [83] using a different RG 
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based techniqu^^ [831 [85] . The implications of this fact are that when we consider 
a one-parameter family of QFTs (i.e. a magnetic system at different temperatures) 
and approach the PT, the critical exponents are determined by the RG properties. 
All the theories on the stable manifold have infinite correlation lengths, and have 
long-wavelength properties similar to the ones of the fixed point theory. 

Furthermore, the existence of two fixed points means that in this theory space 
there exists asymptotically safe theories. In the vicinity of the Gaussian fixed point 
there exists a two-dimensional manifold tangent to the plane {m^. A} consisting of 
theories that in the UV sink into the GFP. We can consider theories that in the IR 
are attracted to the NGFP, so that for them the RG flow is bounded as we change 
the scale, justifying a posteriori our crude treatment of the cutoff A. Of course in 
this particular case this would not be an issue, since this theory is not supposed to 
be fundamental, and we could have a natural UV cutoff A ^ 1/a where a is the 



lattice spacing (nanometers), as well as an IR one of order of the meters. Figure 2.7 
depicts the flow from the GFP into the Wilson fisher one in a truncation including 
up to (j)^ interactions; notice how the initial condition has no six-points interaction, 
which is generated as an effective vertex. 

As mentioned, the FRGE technique can be extended to fermions, gauge theo- 
ries and curved backgrounds, with some technical difficulties A significant appli- 
cation of these ideas is to establish whether Gravity may be asymptotically safe as a 
QFT, due to a UV attractive NGFP, in a scenario similar to the one described here, 
involving however many subtleties [HI [771 CHI [Ml [60l UHl [53] (see [55] for a review) . 

Let us conclude emphasizing that, to us, the most interesting feature of FRGE 
is the possibility of extracting predictions for a class of systems with relatively few 
computational difficulties and, what is most important, knowing only the general fea- 
tures of the systems under consideration, regardless of the details of its microscopic 
dynamics. In this sense, the RG approach explains our success in quantitatively 
describing Nature despite our ignorance of its most fundamental laws. 

Exercise 2.8 (The Ising universality class). Show that 
(2.112) H = -IJ2 J(^' - S{r) S{r') + {S{rf - if , 



r, r' 



that is an Ising-like Hamiltonian, can he written in the continuous limit as 
(2.113) H = j d^x(^-^Ja^R^ S{x)AS{x)-{2K + J)S{xY + KS{xy 



^^The interested reader will find nice pedagogical expositions in [l4l [79l [90 ] , 
^^To mention two: respecting gauge invariance, that can be done by the background field method 
[T], and evaluating traces in curved backgrounds, that requires heat-kernel techniques |651 126] . 
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where a is the lattice spacing and J = J2r '^(^) ' -^^ ~ J2r •^i''^) ■ 

Exercise 2.9 (The local potential approximation). An efficient way of find- 
ing V is, rather than resorting to a truncation, to make an ansatz of the form 

(2.114) TkM= j d''x{-Zkv/\^ + ZkVk{^^)) . 

Show that in this case the fixed point condition is the differential equation for the 
dimensionless potential V{ip) (with rj = 0) 

(2.115) - DVk + (D -2)Vl + - — rju^^ TT ^ ^ = 0. 

^ ' I ; 2(47r)^/2r[2 + Z}/2]l + 2V;' + V^," 

Notice that this can be extended beyond the local potential approximation |2l 
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